Dynamical Tunneling in Mixed Systems 
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We study quantum-mechanical tunneling in mixed dynamical systems between symmetry-related 
phase space tori separated by a chaotic layer. Considering e.g. the annular billiard we decompose 
tunneling-related energy splittings and shifts into sums over paths in phase space. We show that 
tunneling transport is dominated by chaos-assisted paths that tunnel into and out of the chaotic layer 
via the "beach" regions sandwiched between the regular islands and the chaotic sea. Level splittings 
are shown to fluctuate on two scales as functions of energy or an external parameter: they display 
a dense sequence of peaks due to resonances with states supported by the chaotic sea, overlaid on 
top of slow modulations arising from resonances with states supported by the "beaches" . We obtain 
analytic expressions which enable us to assess the relative importance of tunneling amplitudes into 
the chaotic sea vs. its internal transport properties. Finally, we average over the statistics of the 
chaotic region, and derive the asymptotic tail of the splitting distribution function under rather 
general assumptions concerning the fluctuation properties of chaotic states. 



PACS numbers: 05.45. -Hb, 03.65. Sq 

I. INTRODUCTION 

A detailed understanding of how the coexistence of 
classically regular and chaotic phase space areas is re- 
flected in the corresponding quantum dynamics poses 
one of the challenging problems in the field of "quan- 
tum chaos" . Even though semiclassical theories exist 
for the two limiting cases of fully integrable or fully 
chaotic classical dynamics the quantum mechanical 
properties of systems with "mixed" classical dynamics 
have up to date not been amenable to a semiclassical 
formulation. The quest for such a theory is importuned 
by the fact that mixed systems comprise the majority of 
dynamical systems found in nature. 

Out of the wealth of phenomena reported in mixed 
systems, a particularly interesting one is genuinely 
quantum-mechanical in nature: tunneling. A situation 
that has received much attention is the one in which 
tunneling takes place between distinct, but symmetry- 
related regular phase space regions separated by a chaotic 
layer. Interest surged when it was discovered that energy 
splittings can increase dramatically with chaoticity of the 
intervening chaotic layer • This was attributed to a 

suggested mechanism of chaos-assisted tunneling |^,0j8| 
in which tunneling takes place not in a single tunneling 
transition, but in a multi-step process containing tun- 
neling transitions between regular tori and the chaotic 
region, as well as chaotic diffusion inside the chaotic sea. 
Since a large part of the phase space distance is thus tra- 
versed via classically allowed transitions, indirect paths 
can be expected to carry considerably more tunneling 
flux than direct ones. 

Additional evidence was given by the observation that, 
apart from an overall enhancement, the tunneling split- 
tings vary rapidly over many orders of magnitude as a 
function of energy, Planck's constant Ti or other system 



parameters. This was attributed to the occurrence of 
avoided crossings between regular doublets and chaotic 
states which made it possible to further establish 
chaos-assisted tunneling by studying its effect on statisti- 
cal properties such as the splitting distribution function. 
Comparison with predictions of appropriate random ma- 
trix models showed very good agreement 0,^,^- How- 
ever, the lack of a semiclassical description of the tun- 
neling processes remained as a gap between the quantum 
and the classical picture, and — more importantly — the 
size of the tunneling amplitudes was unknown in the sys- 
tems under study, which made a direct and quantitative 
treatment of the phenomenon impossible. 

Both of these problems were addressed in an earlier 
publication by the authors of this work in which a 
semiclassical analysis of tunneling processes in the annu- 
lar billiard was performed, and a formula for the contri- 
bution of chaos-assisted paths to the energy splitting was 
derived. Here, we give a detailed account of our find- 
ings. Particular emphasis will lie on the description of 
how the tunneling rate is affected by phase-space struc- 
tures within the chaotic region, namely the existence of 
an intermediate "beach" region sandwiched between clas- 
sically regular islands and the chaotic sea. 

The structure of this Paper is as follows. In Section ||, 
we review the basic ideas underlying this work. Also, we 
introduce the model system under consideration, the an- 
nular billiard. In Section [II, we introduce the method of 



our analysis, the scattering approach to the quantization 
of closed systems, and explicitly construct the "scatter- 
ing matrix" S for the annular billiard. We then show in 
Section IV how the scattering matrix approach can, un- 
der rather general assumptions, be implemented to the 
study of tunneling in phase space. We explain how S 
can be approximated by a five-block matrix model with 
different blocks representing regular dynamics on either 
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of the islands, beach motion close to each island, and 
chaotic dynamics in the center of the chaotic sea. We 
derive formulas for ^-matrix eigenphase shifts and split- 
tings in terms of paths passing through different combi- 
nations of blocks, laying emphasis on the effects arising 
from the inclusion of the beach blocks. Additionally, we 
track how tunneling flux spreads in phase space and give 
a detailed discussion of the interplay of tunneling proba- 
bilities into, and transport properties within the chaotic 
layer. Finally, in Section ^ we calculate statistical quan- 
tities — such as the splitting distribution function and 
median values for the splitting — by averaging over the 
properties of the chaotic block. We conclude with a dis- 
cussion. 



states can unambiguously be classified as "regular" and 
"chaotic" (for a review, see [^). Regular states are sup- 
ported by classical tori obeying EBK quantization rules 
whereas chaotic states are associated with chaotic 
phase space regions (or subsets of it Q ) . The structure 
of chaotic states is to date not fully understood and is 
presently the subject of intensive research. Classification 
of states as regular and chaotic can become problematic 
at intermediate energy (or h), since EBK-like quantiza- 
tion rules can apply also to states residing on chaotic 
phase space regions lying in close proximity to the regular 
island 16 |]. Loosely speaking, the regularity of clas- 
sical islands can quantum-mechanically extend into the 
chaotic sea, and states of an intermediate nature emerge. 



II. DYNAMICAL TUNNELING 

A. Classical and Quantum Mechanics of Mixed 
Systems 

1. Correspondence of Wavefunctions with Classical 
Structures 

Systems with classically mixed dynamics display both 
regular and chaotic behavior, depending on the starting 
conditions of the trajectory considered. The structure 
of phase space can conveniently be probed by use of a 
Poincare surface of section (PS OS) a phase space 
cut r giving rise to the Poincare map 

iQ,Ph ^ {Q,P)f , iQ,Ph/f = iQ{x,/f),P{x.,/f)) , 

where {Q,P) is a set of canonically conjugate variables, 
and Xi/f G r are connected by the system dynamics. 
If one starts out with highly localized distributions and 
plots iterates of the Poincare cell C-p — {{Q{x), P{x)) : 
X £ r}, then chaotic areas show up as areas which quickly 
become more or less uniformly covered, while regular mo- 
tion remains confined to lower-dimensional manifolds on 
C-p- In a mixed system, both types of structures appear, 
and one arrives at plots of Cp such as the one presented 
in Fig. ^ (see below). 

In order to associate a system's quantum eigenstates 
tp with classical features — such as chaotic regions or 
regular tori — one often uses the Wigner transforma- 
tion [|l2| of the projector \'ip){'ip\. By smoothing over 
minimal-uncertainty wavepackets one obtains the Husimi 
distribution that defines a real, non-negative proba- 
bility density in phase space. We will tacitly invoke the 
Wigner/Husimi concept in the sequel when referring to 
the correspondence of quantum states with phase space 
structures. 

Until now, no general theory for the quantization of 
mixed systems is at hand. However, the understanding 
has emerged that, in the semiclassical limit, quantum 



2. Effect of Phase Space Symmetries 

To discuss the effect of phase space symmetries on the 
structure of quantum states, we consider a system with a 
discrete two-fold phase space symmetry T. We suppose 
that there are two disjoint phase space objects Ai and 
A2 , each of which is invariant under the classical dynam- 
ics, mapped onto another by the symmetry operation, 
A2 ~ TAi- We also suppose that, in the semiclassical 
limit, each of the A1/2 supports a set of states primarily 
localized on it. Let us, for the sake of definiteness, con- 
sider the case when the Ai /2 are EBK-quantized tori. On 

each of the tori, one can construct quasi-modes ipr^^ {q) 
and tpr"^^ [q] — i/'r ('?"q) that obey the Schrodinger equa- 
tion to any order of % The corresponding EBK en- 
ergy eigenvalues Er are then degenerate to any order in 
Ti. However, exact quantum states are constrained to be 
symmetric or antisymmetric under T, 

V'±(g)«-^(^W(9)±^(^Hg)) , (1) 

and the energy degeneracy is lifted by tunneling processes 
by an amount 5Er, giving rise to tunneling oscillations 
with period 2TTh/SEr. 

The best-known example of quantum-mechanical tun- 
neling oscillations is the one-dimensional symmetric dou- 
ble quantum well, where the phase space symmetry 
T{x,p) ~ {—x,—p) connects regular tori in each of the 
wells (for a careful discussion along the above line of argu- 
ment, see [^). In systems of more than two-dimensional 
phase space, symmetries can give rise to more intricate 
situations. The tori Ai and A2 must not necessarily be 
separated by an energy barrier in configuration space, but 
the transition from Ai to A2 can also be forbidden by a 
dynamical law. In this case, there is a dynamical variable 
other than energy that is conserved by classical dynam- 
ics, but violated by quantum dynamics [ p^|Jl9| ]. The case 
of quantum doublets connected by tunneling processes of 
this type was first reported by Davis and Heller pOl who 
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also coined the term dynamical tunneling. A particularly 
clear example of dynamical tunneling will be presented 
in Section [I C in the discussion of the annular billiard. 



As in the case of energy barrier tunneling, splittings due 
to dynamical tunneling can be expected to be very small, 
since classical transport from Ai to Ai is forbidden. 

Note that the formation of doublets is determined by 
the phase space topology of the supporting region, not its 
regularity or chaoticity. The occurrence of doublets has 
also been observed in situations, in which the localizing 
mechanism was due to dynamical localization p6|,p|,pT| 
or scarring Q . Conversely, a phase space structure B ~ 
TB mapped onto itself supports states that do not form 
doublets, regardless of its dynamical nature. 



B. Chaos- Assisted Tunneling 

Apart from the possibility of dynamical tunneling, tun- 
neling processes in systems of more than one degree of 
freedom can have an additional aspect of interest: the ap- 
pearance of chaos in the region of phase space traversed 
by the tunneling flux. As an early paradigm of such a 
system, Lin and Ballentine H) proposed the periodically 
driven double well potential, where chaoticity can grad- 
ually be introduced by increasing the driving strength. 
Lin and Ballentine performed a numerical study of tun- 
neling oscillations between states associated with regular 
tori corresponding to classical motion confined to either 
bottom of the well. They observed that, as the separat- 
ing phase space layer grows more chaotic with increasing 
driving strength, tunneling rates are enhanced by orders 
of magnitude over the rate in the undriven system (the 
integrable 1-d double well). In a later study of the same 
system, Utermann et al. ^ established that the tunnel- 
ing rate of a wave packet initially localized on one regu- 
lar island is determined not by the wave packet's overlap 
with the other island, but by its overlap with the chaotic 
sea, pointing at a role for classically chaotic diffusion as 
a mediator of quantum tunneling flux. 

In a parallel and simultaneous line, Bohigas, Tomso- 
vic and UUmo advocated the interpretation that the en- 
hancement of tunneling was a case of resonant tunneling 
due to the occurrence of avoided crossings of the tunnel- 
ing doublet's eigenenergies with the eigenenergy of a state 
residing on the intervening phase space layer Q . For ob- 
vious reasons, the phenomenon was named chaos- assisted 
tunneling. The interpretation of tunneling enhancement 
in terms of a three-level process was derived from the ob- 
servation that splittings of regular doublets are rapidly 
fluctuating quantities as functions of parameters such as 
energy, Planck's constant ?i, or other model parameters 
— much in contrast to the smooth, at most oscillatory, 
dependence of tunneling rates on h, say, in the case of 
energy barrier tunneling [E3[. 



Since a semiclassical description of tunneling matrix 
elements was lacking, Bohigas and coworkers focussed 
on the statistical fingerprints of chaos-assisted tunneling, 
with emphasis on the consequences of resonance denom- 
inators on the splitting distribution. To this end, the 
interaction of regular doublets with chaotic states was 
formulated in terms of a block matrix model , in which 
properties of states residing on the chaotic sea were ap- 
proximated by use of random matrix theory p4| . This 
model was subsequently refined by Tomsovic and UUmo 
Q to take into account the effect of additional time scales 
in the chaotic dynamics that can appear when residual 
phase space structures, such as cantori, are present in the 
chaotic sea acting as imperfect transport barriers. Pre- 
dictions made using these block-matrix models showed 
good agreement with numerically calculated splitting dis- 
tributions. Along these lines, UUmo and Leyvraz |^ were 
also able to derive analytic expressions for the splitting 
distributions in the case of structure-less chaotic dynam- 
ics, as well as for a structured chaotic sea. Again, the- 
oretical predictions showed good agreement with exact 
numerical data. 



C. The Annular Billiard 

We now introduce the specific system under considera- 
tion in this work, the annular billiard. It was proposed by 
Bohigas et al. and consists of the area trapped between 
two non-concentric circles of radii R and a < R centered 
at {x, y) coordinates O = (0, 0) and O' = (—(5, 0), respec- 
tively. We consider the case oi 5 < a and set i? = 1, 
unless otherwise stated. Note that the billiard is sym- 
metric under reUections at the a;- axis. 




FIG. 1. Parameterization of classical trajectories. Note 
that coordinates ["i,L) are not the BirkhofT coordinates {ip,L) 
usually employed in billiards. The two coordinate sets are re- 
lated by (;9 = 7 -I- Q — TT, where a = arcsin(L/i?). 



1. Classical Dynamics 

Classical motion in a billiard is given as free flight be- 
tween specular reflections at the boundaries. We select 
the PSOS r as a circle of radius r concentric with the 
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outer circle and choose r to be infinitesimally smaller 
than one. Upon in-bound passage through F — or, equiv- 
alently, after reflection from the outer circle — we record 
the trajectory's coordinates {"f^L), where 7 denotes the 
angle of the velocity vector with the x-axis and L — sin a 
is the classical impact parameter with respect to O (sec 
Fig. {'jjL) are canonically conjugate with respect to 
the reduced action p9| ] , and the Poincare cell is given by 
C-p — [0,27r]x] — The billiard's mirror symmetry 

y i—s- — y translates into an invariance of C-p under the 
mapping (7, L) 1-^ (27r - 7, -L). 



correspond to rays along the symmetry axis on the left 
hand side and on the right hand side of the inner circle, 
respectively. 

It will turn out to be of great importance that there 
is a region of chaotic, but relatively stable motion sur- 
rounding each regular island. In the strip of |L| < a + (5 
this stability is easy to understand, as trajectories typi- 
cally encircle the inner disc many times until a hit occurs, 
and at each hit the change in impact parameter is small. 
The "stickiness" of this beach region is increased by the 
presence of regular island chains and of cantori that are 
the remains of broken WG tori. 




FIG. 2. Poincare plot of classical motion in the annular 
billiard at a = 0.4 and 5 = 0.2. 

In Fig. 1^ we present a Poincare plot oiC-p at parameter 
values a — 0.4 and S = 0.2. Rays of |L| > a + S do not 
hit the inner circle, but forever encircle the inner disc at 
constant L filling horizontal strips (or subsets of horizon- 
tal strips) in the Poincare plot. Each of these whispering 
gallery (WG) tori, specified by its impact parameter L, 
is associated with a partner torus — L by the mirror sym- 
metry. These WG tori will be the tunneling tori under 
consideration in this work. 

Rays of intermediate impact parameter \L\ < a + 6 will 
eventually hit the inner circle, and since angular momen- 
tum is then not preserved, motion is no longer integrable. 
This can give rise to the whole range of phenomena asso- 
ciated with non-integrable systems of mixed phase space: 
regular islands and island chains, chaotic regions, partial 
transport barriers (cantori) and the like. The structure 
of the phase space layer \L\ < a + 6 is primarily organized 
by two fixed points of the Poincare map: (i) an unstable 
fixed point at (7,^) = (0,0) and its stable and unstable 
manifolds, along which a chaotic region spreads out, and 
(ii) a stable orbit (7, L) = (tt, 0) at the center of a reg- 
ular island of "libration" trajectories. The fixed points 



2. Quantum Mechanics 

Quantum mechanics of the annular billiard with 
Dirichlet boundary conditions is given by the Hclmholtz 
equation 

(A + fc2)^/.(q) =0 

and the requirement of vanishing wavefunction on the 
two circles. The wave number k is related to energy by 
E — h'^k^/2m. (We note that there exists an analogy 
between quantum billiards and quasi-two-dimensional 
microwave resonator which has proven instrumental in 
many experimental realizations of billiard systems psf.) 

We give here only a qualitative picture of the quantum 
states, deferring a full solution to Section IIIB. It is most 
appropriate to decompose the wavefunction into angular 
momentum components by writing 



(2) 



n— — oc 



where (r, ip) are polar coordinates with respect to O, and 
^n''^\^) denote Hankcl functions of the first and second 
kind, respectively, of order n. We recall that angular 
momentum quantum numbers n are in the semiclassical 
limit related to classical impact parameters L = n/k. 

To understand the nature of quantum states supported 
by the annular billiard, it is instructive to first consider 
the concentric billiard and then to "turn on" the eccen- 
tricity S. If (5 = 0, then angular momentum is conserved, 
and states are paired in energetically degenerate doublets 
composed of angular momentum components n and —71. 
In the eccentric system {S ^ 0), the degeneracy is lifted 
by the breaking of rotational invariance. However, an- 
gular momentum doublets are affected to different de- 
grees — depending on the size of n relative to k{a + 6). 
The symmetry breaking has large effect on doublets of 
small angular momentum \n\ < k{a -\- 5) corresponding 
to classical motion that can hit the inner circle. For low- 
n doublets, the doublet pairing disappears quickly with 
increasing 6 and "chaotic" eigenstates appear that spread 
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out in angular momentum components roughly between 
—k{a+5) and k{a+6). High-angular momentum doublets 
with |n| > k{a + S) are affected only little by the symme- 
try breaking. The doublet-pairing persists, and energy 
degeneracy is only slightly lifted. States are primarily 
composed of symmetric and antisymmetric combinations 
of n and —n angular momentum components, 

|«(±))«-L(|n)±|-n)) . (3) 

Note that each of the quasi-modes | ± n) corresponds to 
classical motion on the WG torus ±L. 

We present plots of a "regular" doublet quantized at 
/c ~ 55 and a "chaotic" state at A: « 60 in Fig. ^ and 
compare them to classical trajectories with starting con- 
ditions in the chaotic sea and on WG tori, respectively. 
The correspondence between quantum states and the na- 
ture of classical dynamics is clearly visible. 




FIG. 3. Comparison of classical dynamics and quantum 
states of the annular billiard, (a-c) Classical motion: (a,c) 
regular trajectories of ±L, \L\ > a + 5, (b) chaotic trajectory 
of \L\ < a + S. (a-c) Eigenmodes: (a,c) doublet of "regular" 
eigenstates at ~ 54.434, and (b) "chaotic" eigenstate at 
k ^ 60.252. 



3. Tunneling Between Whispering Gallery Tori 

Let us discuss the high-angular momentum doublets 
in more detail. As explained above, the energy splitting 
6 En between |a(+)) and \a^~'>) gives rise to tunneling 
oscillations between quasi-modes | ± n) associated with 
WG tori ±L = ±n/k {L > a + S). A quantum particle 
prepared in state \n) will therefore change its sense of 
rotation from counter-clockwise to clockwise and back to 
counter-clockwise with period iirfi/SEn. Note that this 
tunneling process serves as a particularly clear example 
of dynamical tunneling. It occurs in phase space rather 
than configuration space, as the corresponding tori are 
identical in configuration space. Also, the tunneling pro- 
cess does not pass under a potential barrier in configu- 
ration space. In fact, energy does not play any role in 
the tunneling, as energy is related only to the absolute 
value of the momentum vector and not to its direction. 
Rather, the tunneling process violates the dynamical law 
of classical angular momentum conservation for rays of 
large impact parameter. 

The concept of chaos-assisted tunneling can be nicely 
visualized for the case of tunneling between WG modes 
in the annular billiard. In fact, the annular billiard was 
proposed as a paradigm for chaos-assisted tunneling in 
Ref. 01 . In chaos-assisted processes, tunneling tori ±L 
are connected not by direct transitions between n i— )■ — rt, 
but by multi-step transitions n i— s- t . . . — I' i— > —n. 
A particle tunnels from n to some i ^ k{a + 6), tra- 
verses the chaotic phase space layer by classically allowed 
transitions to reach the opposite side of the chaotic sea 
—i' ^ —k{a + 6), and finally tunnels from there to —n. 
To establish this notion, Bohigas et al. [Q checked the 
behavior of splittings as the eccentricity is changed to 
make the intervening phase space layer more chaotic. In 
a numerical study, they compared the splittings of reg- 
ular doublets with the rate of classical transport across 
the chaotic layer. The findings showed that the splittings 
increase dramatically over many orders of magnitude as 
chaotic transport becomes quicker. 

However, without a quantitative — possibly semiclas- 
sical — theory, it is impossible to separately analyze the 
importance of tunneling amplitudes and transport prop- 
erties of the intervening chaotic layer. Usually, the pa- 
rameter governing symmetry breaking in a mixed system 
changes both tunneling amplitudes from/into the regular 
torus and chaoticity in the intermediate layer. In or- 
der to separate the relative importance of these effects, 
a quantitative understanding of the tunneling processes 
must be obtained. Such a quantitative description of 
chaos-assisted tunneling was given in Ref. [|lO| and will 
be developed in full detail in the sequel. 
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III. QUANTIZATION BY SCATTERING 

A. General Description of the Method 

In this work, we will employ an scattering approach 
to quantization |2^j27[] which, in essence, is constructed 
as the quantum-mechanical analogue of the classical 
Poincare surface of section method. For the sake of self- 
containedness, we give a brief review of the scattering 
method. 

Let us consider the case of a billiard G and introduce 
a Poincare cut F in configuration space, thereby dividing 
Q into two parts, ^7+ and Q^. We suppose that F can be 
chosen along a coordinate axis (the (72-axis, say) and that 
the wave problem is separable on an infinitesimal strip 
around F. At a given energy, one chooses a complete set 
of functions (j)n (92) along F and writes the wavefunction 
on an infinitesimal strip around F as 



the (fin''^'^ (j)"n' so that in the region of classically al- 
lowed motion, the contribution is exponen- 
tially small for all but a finite number of indices (the 
so-called "open channels"). This allows the truncation 
of S to finite dimension, say \n\ < A, with an error that 
is exponentially small. Both scattering matrices S± can 
be constructed in a representation such that they are 
unitary in the space of open modes and, if the system is 
time-reversal invariant, symmetric. S on the other hand 
is unitary, but not necessarily symmetric. In spite of this, 
we will in this paper also refer to S* as a scattering matrix. 

It is clear by construction that S is the quantum- 
mechanical analogue of the Poincare mapping ||2^. Its 
iV-th iterate constitutes a time-domain-like propa- 
gator. Note that the iteration count N of the Poincare 
map does not correspond to a stroboscopic discretization 
of time, but rather to a, fictitious discrete time, since gen- 
erally the time elapsed between passages of F can vary. 
If, however, we select our surface of section properly, then 
return times will not vary by too much, and S will not be 
too different from the genuine time-domain propagator. 



and (pn' •* are functions that, in the semiclassi- 

cal limit, correspond to waves traversing F in positive 

and negative gi-direction, respectively. We assume that 
f 2) 

the set (j)n ((72) is chosen such that quantum numbers n 

(2) 

correspond to values kn of longitudinal wave number. 

Then, fci^^ and transverse wave numbers kn"^ are related 

by E„ ^ a^[(fci^^)2 + {kl?)'^]/2m. Note that orthogonal- 

f 2) 

ity of the modes on F is ensured by the choice of the ' . 
Each of the domains Q± constitutes a scattering system 
that scatters waves (pn'^^ into waves (j3^n'~^ and 
vice versa. Associated with these scattering systems Q± 
are scattering matrices S±{E) that relate the coefficient 
vectors 



/3 = S-{E)a. 



and 



a = S+{E)I3 



(5) 



The quantization condition is equivalent to the require- 
ment of single- valuedness of the wave function on F, and 
so to the equivalence of the two scattering conditions. 
The system supports an eigenstate whenever the prod- 
uct matrix 

S{E)^S^{E)S+{E) 

has an eigenvalue of unity, hence the quantization condi- 
tion reads 



det{S{E) - 1) = . 



(6) 



At a quantized energy, the wave function can be recon- 
structed from the corresponding eigenvector a of S' via 
Eqs. In principle, S is an infinite-dimensional ma- 

trix. However, in many cases of interest one can choose 



B. Scattering Matrix of the Annular Billiard 

It is fairly straightforward to implement the scatter- 
ing approach to the case of the annular billiard. As 
discussed in Section IIC, we choose F as a circle of ra- 
dius r, where a -\- 5 < r \. Since classical impact 
parameter is conserved by motion on the WG tori, we 

(2) 

choose q2 = (fi and (f)n — exp{in(p) on F. Waves travers- 
ing F are given by outgoing and ingoing cylinder waves 
4^'+) = i"H,W(fcr), (t>i^'~'> = i"H,(2)(fcr), and we ob- 
tain the decomposition in Eq. (0), 



V'(r,(p)= i" Unili^\kr) + /3r,lli^\kr) 



imp 



(7) 



n— — 00 



In the present example, outgoing waves are scattered 
to ingoing waves by the interior of the outer circle — 
giving rise to the scattering condition a = S' [k) (3 — 
and ingoing waves are reflected off the exterior of the 
inner circle, which leads to the relation /3 = 5' (fc) a . In 
order to derive the explicit formulas for S '^'^^ and S , we 
note that the Dirichlet boundary condition ijj{R, (p) = 
on the outer circle leads to 



(8) 



S (k) is diagonal, in accordance with angular momen- 
tum conservation in scattering events off the outer circle. 

S (k) is derived by performing a coordinate change 
to the primed coordinates defined with respect to O' = 

0). In this representation the scattering matrix is 
S"(i)(fc) = -H^ )(fca)/H^i'(fca)5„,™. Transformation 
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back to unprimed coordinates is done by the addition 
theorem for Bessel functions (see e.g. pof) 

C30 



for (5 < a = r'. J„(.t) = (H^^^(a;) + lli^\x))/2 denotes 
the Bessel function of order n. We arrive at 



For eigenvectors a^^' of S*, 



«=-oo 



Time-reversal invariance imposes the symmetry 
S^^nik) = on S^^Hk) and by virtue 

of the mirror symmetry S'l'^ _„j(fc) — Sn}n{k). (Note 
that the formula given for S^^^{k) differs by the factor 
i ("^'") from a formula given earlier by us and con- 
sequently, also the symmetries are different. This is due 
to a slightly different choice of basis in (^.) The full 
5-matrix is then given as 



0) _ . (i) 



(10) 



where aj € {±1}- Whenever the system supports an 
eigenstate, aj determines the symmetry of the corre- 
sponding wave function with respect to the x-axis. 

In Fig. ^, we show a grayscale plot of |S'„,m| as a func- 
tion of ingoing and outgoing angular momentum for the 
parameter values a = 0.4, 5 = 0.2, k = 100. The overall 
structure of S is governed by classically allowed transi- 
tions: it is mainly diagonal in the region of high angular 
momentum |m| > k(a-\-d), whereas the inner block re- 
flects the dynamics given by the classical deflection func- 
tion. Here, the main amplitude is delimited by two ridges 
that correspond to classical rainbow scattering. 

The tunneling amplitudes relevant to the WG splitting 
are contained in S as non-diagonal entries Sn^m with n > 
k{a + S). Fig. H depicts |S'„,m| for the above parameter 
values and n = 70 (at k{a + 6) = 60). 
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FIG. 4. Grayscale plot of the quantum 5-matrix |5'n,'Ti| 
(arbitrary gray scale). The dashed line indicates the region 
of classical angular momentum mixing |n| < k(a -I- S). 

Using the relation Yi^}f^{x) = (-)" H^^'2)(x) for in- 
teger n, one verifies that the spatial symmetry of the 
annular billiard translates into the S'-matrix symmetry 
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m 

FIG. 5. Absolute values of tunneling matrix elements 
l-^n,™! for n = 70 and m — —100, . . . , 100 (logarithmic scale). 
Inset: first quadrant of S, showing the tunneling ridges. 
(Taken from Ref. @.) 

The tunneling amplitudes are largest 10^^) around 
m = 63 and fall off faster than exponentially away from 
this maximum. As can be seen in the inset of Fig. 
the line of maximal tunneling amplitudes continues the 
line of rainbow ridges into the regime of classically for- 
bidden transitions. Close to these tunneling ridges, and 
in the direction away from the diagonal, one can observe 
oscillations in jiSn^ml that resemble the Airy oscillations 
well-known in diffraction theory | |30|] . It is interesting 
to note that the behavior of the tunneling probabilities 
is not a monotonous function of |n — m|, that is, of the 
phase space distance traversed. A mere distance in phase 
space can therefore not serve to estimate the behavior of 
tunneling probabilities. 
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We finally mention that the tunneling transitions con- 
sidered here can be given a semi-classical interpretation 
in terms of complex rays that interact with the analyt- 
ical continuation of the inner circle to complex config- 
uration space jl^j3l|,^. These rays can either scatter 
off the inner complex circle by a generalization of spec- 
ular reflection, or they can creep along a complexified 
inner circle (of complex radius determined by the poles 
of the internal scattering matrix) by a mechanism simi- 
lar to that proposed by Franz 1 33 1 and Keller . Every 
tunneling pair , n/ is connected by at least one com- 
plex reflection trajectory with real initial (final) impact 
parameter Li = rii/k {Lf = Uf/k) and complex initial 
(final) angle 7^ (7/)- Complex reflection and creeping 
trajectories give rise to contributions to the tunneling 
matrix elements that scale with k like 

■Srofloction ~ exp - fcArcf(«, /) ) , 

5'crccping - gxp ( - kA„{l, f) + k'^^B„{l, /) ) , 

respectively, where A^d, ^cr and B^r are classical prop- 
erties of the associated trajectories that are independent 
of k. The relative importance of the two contributions 
depends on fc, the billiard geometry, and the initial and 
final angular momenta considered. At the present param- 
eter values, the contribution due to refiected rays usually 
dominates the one arising from creeping rays. 

The semiclassical picture provides an intuitive expla- 
nation to the tunneling ridges mentioned above: they 
can be identified as combinations n, m where one of the 
angles 7i,7/ is closest to reality. Also, the nearby oscil- 
lations can be understood as arising from a coalescence 
of two reflection saddle points. 



IV. TREATMENT OF CHAOS-ASSISTED 
TUNNELING 

A. Implementation of the Scattering Approach 

We now discuss how the scattering approach can be im- 
plemented to the treatment of chaos-assisted tunneling. 
Let us suppose that the system under consideration has 
a phase space symmetry of type ((72,^2) '-^ (^ <Z2,— P2), 
and that the g2-axis can be chosen as a PSOS F. We as- 
sume that the wave problem is locally separable around 
F, which renders the problem tractable by the scatter- 
ing approach. Let us now pick two classical objects Ai 
and A2 = TAi and consider quasi-modes 'tpr^^'^\q) sup- 
ported by Ai/2, as discussed in Section ||. We assume 
that motion on Ai and A2 corresponds to conservation 
of p2 and —p2, respectively, and that ±p2 is semiclassi- 
cally related to the quantum number ±n. (Note that in 
a general system, the choice of the PSOS F and a proper 
basis for the scattering matrix can be a very difficult 



task. In this Paper, we will not deal with the problem of 
solving a general scattering problem — except for the an- 
nular billiard — and assume the scattering matrix Sn^m 
to be known.) By symmetry p2 1— > — P2, we find that 
Sn,7i — S-n,-n, and by virtue of the classical conserva- 
tion of p2, Sn,n is almost unimodular. The deviation of 
|'S'n,n| from unity will be due to classically forbidden (i.e. 
tunneling) transitions, and so is expected to be small. 

Let the quantization energies of the doublet be de- 
noted by E^, and let \a'^{E^)) be the eigenvectors corre- 
sponding to the two quantum states. We now make the 
approximation that the properties of these two vectors 
are, to good precision, given by the eigenvector doublet 
\a.^{E)) at one fixed energy E lying between E'^ and 
E^ . The corresponding (generally non-zero) eigenphases 
be denoted by 9^{E). 

Dropping the energy variable E, we decompose the 
eigenphases 9^ in the form 



i^) ± \ 



(11) 



where A9n^'^^ are taken to be real, and = 
— i logs',!, „. The quantities A9n = A.9n^^ 4- i and 
69n can be interpreted as the shift and the splitting, re- 
spectively, of the exact eigenphases due to tunneling pro- 
cesses. These eigenphase quantities are trivially related 
to the energy shift and splitting, as will be explained be- 
low, and it is therefore sufficient to calculate A9n and 
S9n- Note that the 9^ are real by unitarity of S, and 
therefore also ^i"'' -I- A6'„ and S9n are real quantities. We 
also note that 



1±,± 



e'^^" explzTV 



A9n ± - S9„ 



(12) 



for any integer N. 

From the eigenvector doublet \q.^) we can now obtain 
the vector equivalent of quasi-modes 



±n) 



a. ' 



)± I a 



(13) 



It is clear that \n) and | — n) are localized at (or around) 
components ±n, respectively. Using the symmetry of S 
and Iq;^), we write 



(14) 



In order to derive formulas for A9n and 59n that 
relate these quantities to matrix elements of S^, we now 
combine Eqs. ( |l2|) and (|lj). We choose a positive N 
which satisfies 



iV« \A9,,±59,J2\- 



(15) 
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and expand the second exponential in Eq. (y_2|) to first 
order. Considering upper signs in Eq. ( [1^ ) and taking 
imaginary parts we obtain 



Similarly, taking lower signs in Eq. (jl^) gives 
.{e-'^'^'^"'(n|5^|-n)} . 



2 

N 



(16) 



(17) 



It is instructive to rephrase Eq. ( p7| ) for the splitting: 
starting from Eqs. (|l^,|l^) one can also contract the ex- 
ponentials to a sine and arrive at 



sm 



= |(n I 5^ I n 



(18) 



which has the form of tunneling oscillations in time, 
\sm{5Et/2fi)\ — |(n| exp(— i_fft/?i)| — n)\, with iteration 
count N taking the role of time and eigenphase split- 
ting taking the role of energy splitting. By considering 
Eq. ( p^ ) we therefore probe the onset of tunneling oscil- 
lations in the linear regime. 

Use of Eqs. ( p^ , p^ ) for low N will however necessitate 
precise knowledge of the eigenmodes | ± n) (see e.g. ||3^] 
for a recent application of a similar formula for = 1). 
However, to exponential precision, eigenvectors may be 
just as difficult to obtain as the shift or the splitting itself. 
It is therefore important to realize that use of Eqs. ( p^|jl7| ) 
for large N may allow to extract the quantities of interest 
using much less precise eigenvector information. Instead, 
one then uses dynamical information — in the framework 
of time-domain like propagation with S — which will 
eventually allow the interpretation of tunneling processes 
in terms of sums over paths in phase space. To this end, 
let us reformulate Eq. (|l^) by writing 



(19) 



where the are expected to be small and k_ 
according to the symmetry of |q!*). The right-hand side 
of Eq. (|l3) then reads 



2(n I I ±n 



2 [S 



J n,±n 



(20) 



Here, [S ]n,±n denotes a matrix element of the A^-th 



iterate of S, and C 



with 



+ Ei-tX-t [n™,™. • (21) 

rn.m' 

We see that (nlS*^! ± n) can be replaced by (n|S'^| ± 
n) at the price of corrections 0{k) at most. However, 



since the left hand side of (20) is of size sin(A^A6'„/2) ^ 
0(1) (considering the positive sign) for large N, the first 
term on the right hand side must also grow with TV and 
become of 0(1). In particular, for N > k/\A9^^^\, the 
K-corrections can be neglected. The same holds for the 
negative sign in (po|), but then N ^ K/\S9n\ must hold. 
Therefore, whenever either of these lower bounds on N 
holds simultaneously with the upper bound ( p^ we can 
write 



Ae^«^^^m{e--'^^"'[n„.„}, (22) 
.{e--'^^"'[n „,_„}• (23) 



2 

N 



We note that the conditions on N necessary for Eq. ( p^ ) 
are met if |A6'„| > \69n\ which, as we will find later, is 
always the case. However, the conditions for Eq. (|2^ ) 
might not be met simultaneously. In this case, one has 
to expand in (|l^) to obtain 



2 

N 



[s 



Nl 



(24) 



requiring k/|(56'„| <C <C |^0n|~^ — a condition that 
can always be fulfilled (if k <^ 1). However, we will in 
the sequel calculate the splitting by use of Eq. (|2^), as 
this expression constitutes a linear relation between the 
splitting and the different contributions. We therefore 
assume that the use of Eq. (|2^) is justified and only com- 
ment on the use of Eq. (|2|) . (We note that Eq. (|2|) was 
employed in an earlier account of this work [|lO|.) 

Returning to Eqs. (|2^j2^) recall that any matrix ele- 
ment [S^]n,m can be expressed as a sum over paths in 
matrix element space of length N that start at n and 
end at m by writing out the intermediate matrix multi- 
plications. This yields the real part of the shift and the 
splitting as 

N-l 

Ae^«lsm{e-^«i"' ^ n^A.a.,.}, (25) 



TV 



^m < e 



{n — >n} i—1 

N-l 



{e-^«i"' y: n (26) 



The eigenphase shift and the splitting are therefore given 
in terms of paths of length iV that lead from index n 
back to n or to — n, respectively. Note that in order to 
contribute to the shift, the path must leave the index n 
at least once; the trivial path of constant matrix index n 
does not contribute, as exp(— iiV^i"^) [S^]n,7i = 1 is real. 

For the sake of completeness, we also comment on the 
imaginary part of the shift. Since 6'i°'' -I- A6'„ is real. 



(I) _ 



Smj^i"^}. By unitarity of S, one finds 



A0(^) : 

n 



2 5^ |'5'n,m| 



(27) 
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It remains to connect eigenphase shifts and splittings 
to the corresponding energy quantities. At a given en- 
ergy E — which need not necessarily be an eigenen- 
ergy of the system — the eigenphases 9j{E) are dis- 
tributed on the unit circle, with the position of regu- 
lar doublets determined by Eqs. (pi], [25|, p6|) . Regular 
doublets revolve around the unit circle with "velocity" 
given by the energy dependence of 
the On^\E). These quantities can be taken constant on 
the scale of the energy splittings, and eigenenergy split- 
tings are therefore trivially related to eigenphase split- 
ting. Consequently, we will in the sequel consider shifts 
and splittings of eigenp/iases rather than eigenenergies. 
This has the numerical advantage that we can consider 
S{E) at any E without having to worry about quantiza- 
tion. 

In order to extract the shift or the splitting of the dou- 
blet peaked at indices ±n wc therefore have to consider 
paths in 5-matrix index space that lead from n back to n 
or to — n, respectively. Let us briefly focus on the split- 
ting. In principle, there are two types of paths: direct 
paths and chaos-assisted ones @J^,01. Direct paths tun- 
nel directly from n to —n in a single transition over a long 
distance in phase space. Their contribution to the split- 
ting is of the order IS'n^-.nj. Chaos-assisted paths include 
at least two tunneling transitions over relatively small 
phase space distances. They tunnel from n to some index 
£ such that £ lies in the inner block of classically chaotic 
motion, then propagate — via classically allowed tran- 
sitions — to some £' within the inner block and finally 
tunnel from £' to n. Contributions arising from chaos- 
assisted paths are then of the order \Sn.e Si' ^-n\- As we 
have seen in the example of the annular billiard, tunnel- 
ing matrix elements Sn,£ fall off very rapidly (faster than 
exponential) for large phase space distances |n — £|. This 
rapid decay strongly suppresses the contributions from 
direct paths and explains why the combination of two 
tunneling transitions can be much more advantageous. 



B. A Block-Matrix Model 

We now formulate a generalization of the block matrix 
models usually encountered in the treatment of chaos- 
assisted tunneling that takes into account the effect 
of the transition region between classically regular and 
chaotic motion. 

Statistical modeling of chaos-assisted tunneling is usu- 
ally done in terms of a block matrix model of the type 
regular- chaotic-regular in which properties of the chaotic 
block are approximated by random matrix ensembles . 
This three-block approximation, however, discards all in- 
formation about phase space structures inside the chaotic 
sea, such as the inhibition of mixing by broken invari- 
ant tori (cantori) p&. Inhibition of classical transport 



can lead to dynamical localization of states in regions of 
phase space. Chaotic states then do not extend over the 
full chaotic region of phase space any more, but only over 
components of it. This additional structure would not be 
reproduced by the approximation by a single random ma- 
trix block. Block matrix models for chaos-assisted tun- 
neling were amended to the presence of imperfect layers 
in Refs. j8|J|] by introduction of separate, weakly coupling 
blocks for each of the phase space components. A rela- 
tionship between classical flux crossing imperfect trans- 
port barriers and quantum Hamiltonian matrix elements 
was given in Ref. 

We now argue that the treatment of chaos-assisted tun- 
neling in a generic mixed system usually requires a five- 
block model at least. The reason is that classical motion 
in the "beach" regions close to regular domain is rela- 
tively stable, despite of its long-time chaotic behavior. In 
particular, transport in the direction away from the regu- 
lar phase space region can be strongly inhibited. This dy- 
namical stability leads to the formation of quantum beach 
states that have most of their amplitude in the beach re- 
gion and little overlap with the chaotic sea. It has been 
reported on many occasions that beach states have great 
similarity to regular states residing on the adjacent is- 
land and that they follow EBK-like quantization rules 
]l5| , ^ 16,^. Note that the importance of the beach re- 
gion is highlighted in chaos-assisted tunneling processes: 
as tunneling amplitudes Sn,i decay rapidly away from the 
regular island, chaos-assisted paths of largest amplitude 
will typically lead to indices £ and £' such that the cor- 
responding momenta P2{£) and P2{£') he just inside the 
chaotic sea on either side, that is, in the beach regions. 

In order to take account of the special role of the beach 
regions (or "edge" regions, we will use these two expres- 
sions as synonyms), we propose to generalize the usual 
three-block model regular- chaotic-regular to a five-block 
model of the type regular- edge- chaotic- edge-regular. (In 
this work we assume that, apart from the edge layers, 
no further transport-inhibiting structures are present. 
The existence of further transport-inhibiting structures 
inside the chaotic sea will require the addition of further 
blocks.) 

In the sequel, we approximate S' by a flve-block model 
S as depicted in Fig. ^in which each regular region, each 
beach region and the center chaotic region are modeled in 
a separate block, and coupling between different blocks 
is weak. We assume that S has, by a unitary transforma- 
tion, been converted such that all intra-block transitions 
vanish. We use indices n and — n' for the properties of 
the two regular blocks, £ and —£' for the beach blocks and 
7 for chaotic states. For the diagonal elements, we write 
^A,A = exp(i0A), where A € {n, ^, 7, — — n'}. Note the 
symmetries Sn,n = S-n.-n and S-i^-i = Si^i, and there- 
fore 6^n = Sm 0^1 = Oi- Inter-block coupling elements 
are denoted by Sn.e, Sn,-y, Se,ry and so on. It is natural to 
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assume that the tunneHng elements Sn,e between regular 
tori and the beach region will be much smaller than the 
transition amplitudes Sg^^ between beach and the center 
block. 




chaotic 7 

beach -t 
regular -n 



FIG. 6. The modified block-matrix model: 
blocks participating in the recer contributions. 



structure of 



coupling matrix elements Sn^j, Sg^^ in terms of the orig- 
inal S'-matrix elements — or, even more, as a semiclas- 
sical expression — is presently not possible, as too little 
is known about the precise nature of the quantum local- 
ization on the beach layer. Also, standard semiclassical 
methods break down in the beach. 

For the brevity of notation, we will drop the tildes on 
S and e. 



C. Extracting the Shift, Splitting and Eigenvector 
Structure 

We now use the block matrix model and formulas 
( ^ , p6|) to extract approximations for the shift and the 
splitting of eigenphase doublets. Also, we comment on 
how to extract eigenvector structure and to calculate 
eigenphase properties from it. 



We now explain how the properties of S can be ex- 
tracted from the original matrix S. Note that a block 
that is almost diagonal S will be changed only little in 
the transformation to S . The outermost regular blocks 
are diagonal in all orders of 1/fc (or K), and matrix ele- 
ments of the transformed matrix will differ from the orig- 
inal ones only by exponentially small corrections. We can 
therefore approximate S'„^„' ~ Sn,n^n,n'- The argument 
also holds for the edge region, since the choice of basis for 
the regular blocks will also be good in the edge blocks, 
and we can approximate Se^ei ~ Se^i Se^e' and in the non- 
diagonal block Sn^i ~ Sn,i- Consequently, the choice of 
the border index between regular and edge blocks does 
not affect the results to the present approximation. 

We model quantum dynamics within the inner block 
by a superposition of two Gaussian ensembles p8| , ^ , p8| 
— in our case two Circular Orthogonal Ensembles (COE) 
— that approximate the sets of chaotic states with 
even and odd symmetry, respectively. For a given reg- 
ular mode n coupling matrix elements Sn.-y are chosen 
as Caussian distributed independent random variables of 
variance determined by cr^ ^ = |<Sn,cP/(2^cOE + 1) with 



^COE 

E 

S=-^COE 



(28) 



where we have taken the chaotic block to extend from 
^COE to —icoE- Coupling matrix elements Si^j are de- 
fined analogously. 



15, 



^COE 

E 

9=-^COE 



(29) 



The choice of f coe contains some uncertainty that will 
affect the results as an overall factor in the effective cou- 
pling elements. Note that an approximation for single 



1. Eigenphase Splitting 

When using Eq. (|2^) and the block matrix representa- 
tion of S to calculate the splitting, we have to perform 
the sum over paths 



V. 



N 



E 



n ^ 



{n ^—n} i=l 



of length TV ^ K,/\59n\. As N is large (and under the 
assumption that none of the internal diagonal elements 
S\i,\i is equal to Sn,n) it is sufficient to collect contri- 
butions that are of order N . We consider three families 
of paths: paths of type regular-regular (rr) that contain 
one single jump from n to —n, paths of type regular- 
chaotic- regular (rcr) that pass through the center block, 
and paths regular- edge- chaotic- edge- regular (recer) that 
pass through all five blocks. Adding up the different con- 
tributions gives the eigenphase splitting as a sum 



(30) 



Within each of these families, one has to sum over indices 
of staying times of the path at each diagonal element. We 
show in Appendix A how this can be done and quote the 
result: The sum over paths with M -\- \ jumps passing 



through M intermediate blocks 5^*^'*^-', . . . , S'^ 



iM,iM) 



given as 



-,N (ii,...,iM) 



E 

Ai,...,A, 



M 



5a. 



(31) 



where the sums over the run over all indices of the 
corresponding blocks S'^*'"''"', and \m^\ = —n. Correc- 
tions to (31) are of lower order in N or higher order in 
transition amplitudes. The phase denominators 
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j(0) 



= d 



arise from the summation over staying times Ti at the 
different blocks. Within a given family, each path con- 
tributes a factor eyip[iTi{6 \^ — 6*,!"^)], and geometrical 
summation over Ti results in the denominators listed. 
Note that only phase differences appear that combine 
the outermost phase On^ with one of the phases of 
the inner blocks. Contributions containing other phase 
denominators decay exponentially in N . Also, we have 
only taken into account paths that pass through each of 
the inner blocks once. This "never look back" approx- 
imation is justified since paths containing loops are of 
higher order in the inter-block transition elements. For 
a treatment of loops in index space, see Appendix A. 

Let us now discuss the contributions of the different 
families in turn. For paths of type regular-regular, we 
apply Eq. dp) for M = and find 



1} Sn.—n 



(32a) 



Chaos-assisted paths of type regular- chaotic-regular visit 
the center block. By application of Eq. (^Tj) for M = 1, 
we find that the rcr-contribution to the splitting is 



1 (rcr) 



2 3m < e 
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(32b) 



Finally, the paths regular-edge-chaotic-edge-regular pass 
through three intermediate blocks (Af = 3), hence 



Sn,e J Sj S. 



f dn.i dji '^ dji 



p, J 



(32c) 



One sees that all tunneling rates mediated by internal 
blocks can be strongly enhanced by two effects. 

(1) Combinations of tunneling matrix elements may 
become progressively more advantageous as more steps 
are allowed. 

(2) Coherent summation over staying times results in 
phase denominators. Avoided crossings of these phases 
turn the splittings into a rapidly fluctuating quantity 
with respect to small changes in energy, say, or an ex- 
ternal parameter of the system. The phase denominators 
also lead to an overall increase in the tunneling rate since, 
at a given wave number fc, there are of order k internal 
states available. This means that typically, there will be 
one phase denominator of size d~\^ ~ fc/27r at least. 

Both effects, (1) and (2), can also enhance the re- 
cer contributions with respect to the rcr ones. We can 
therefore expect the recer contributions to dominate the 
tunneling rate. For a given system, their relative im- 
portance may vary, depending on the size of \Sn,'y\ and 
\Sn,eSe^'j/dn,e\- 



2. Eigenphase Shift 

Paths that contribute to the real part of the shift lead 
from n back to n and have to leave this index at least 
once. To do so, they can either tunnel to the center block 
or to the edge block, which leads to a decomposition of 
contributions into 

The summation over these paths can be done by the same 
procedure used above, and one finds that the contribu- 
tions to the real part of the shift are 



^0 (rcr) 



\^ " ^ d^J / ' 



ym < e 



ym < e 



^ d^i J 



E 



'n,7 

Sn,£ Si^j S^J' Sf/ 



(33a) 
(33b) 
^} .(33c) 



Due to the rapid decay of tunneling matrix elements we 

(rcr) I 



expect that |Afti'"''')| » 



^1 > \A0rr'''\. Fur- 
thermore from Eqs. (H^^i), | A^i'''^"'^'^) | ~ {SelT'""'''^ \. 
Consequently, the shift will typically be much larger than 
the splitting. 



3. Eigenvectors 

We can also use the block matrix model to approximate 
regular eigenvector doublets a* . We set = 1/ \/2 and 
Q!^„ = ±l/\/2, and solve 

to leading order in the coupling matrix elements between 
neighboring blocks (neglecting all other coupling matrix 
elements). We find that the components of in the 
beach regions are 



1 Sn,t 

V2 dni 



(34a) 



and a^g — ^or^ ■ Components in the center block are 
given by 



= ^7^ 



\/2 £' S-f'-r 



(34b) 



if the symmetry cr^ G {il} of the block-diagonalizing 
vector 1 7) is the same as that of a*, and = 
otherwise. In both Eqs. ( p4[ ) the relative error is 

Oi\SnAVdn,iASi,^\yd n,7 J ■ 
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This approximation for the eigenvectors can be used 
to relate eigenphase properties to the size of eigenvector 
components in the different blocks. Upon comparison of 
and (34a), we see that the dominant contri- 



(35) 



Eqs. (g3 

bution to the shift can be written as 



where a is either of the or^. Therefore, the eigenphase 
shift is related to the eigenvector's overlap with the beach 



region. Similarly, from Eqs. (B2c) and (34b), 



ecer) ^ cj^^^ . 
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d„,^ [{a+f - [a-f]] . (36) 



Eq. ( ^6| ) relates the presumably dominant contribution 
to the splitting to the eigenvector's overlap with the cen- 
ter block. This explains and quantifies the observation 
of Utcrmann et al. |6| , that regular doublet splittings in 
a mixed system are in close correlation with the states' 
projection onto the chaotic sea. Similarly, Gerwinski and 
Seba |^l[ related tunneling rates between a chaotic phase 
space region and a regular island to the overlap of a 
chaotic scattering state with the regular island. How- 
ever, we see that the ad hoc association 



El 



(37) 



is not complete: in the exact relation (36), each summand 
is weighted by a phase difference. 



4- Comments 

In view of the explicit formulas, we see that the re- 
sults are not significantly changed by approximating the 
regular and edge blocks of S by the corresponding ele- 
ments of the original matrix S. Only in the immediate 
vicinity of resonances between regular and beach eigen- 
phases the approximation Si^i ~ Se^e is not appropriate, 
as it over-estimates the imaginary part of the phase 9i 
and leads to a spuriously broad resonance. (With the 
neglect of A0^^^ we shall not be concerned, because we 
do not aim at an exact reproduction of the peak posi- 
tions.) For most beach states, it is therefore more appro- 
priate to make the somewhat ad hoc approximation of 
Be = a.rg{Sg^g} instead of — ilog5^^^. This approximation 
now leads to an under-estimate of the resonance width, 
but we have checked that in the cases discussed in Sec- 
tion IV D the edge £ are sufficient large {£ > 59) that the 

and \Sn,n — '3iip{i9i)\~^ cannot 



resonances \Sn,n ~ Se,_, 
be distinguished by eye. 



From a methodological point of view, it is worthwhile 
mentioning that the formulas for the shift and the split- 
ting can also be derived from a complementary approach. 
One can expand the characteristic polynomial 

Ps{x) = det(S'- a;) 

of S around Xn = exp{i6}P^) to second order in the ex- 
ternal coupling elements |5'„_Ai |, \SxM,-n\, and then solve 
for its roots x^, Ps{x^) = 0. Upon definition of shift and 



splitting via x 



Ax ± dx/2, one finds formulas 



that, to lowest order in the in ternal coupling elements, 
are identical to Eqs. (32c, 33a). Hence, the "never look 
back" summation over paths corresponds to the lowest 
order of a formal expansion of the eigenvalues, contain- 
ing the coupling elements as small parameters. 

It is important to realize that the different phase de- 
nominators dn,e and dn^-y may fluctuate on different scales 
as functions of the energy or an external parameter. Typ- 
ically, the center block will be much larger than the beach 
blocks, and avoided crossings of 61,1°^ with one of the 9^ 
will occur more often than those with one of the 6i. Also, 
since beach states display EBK-like behavior with actions 
that can be similar to those of the regular states, we can 
expect phase differences 9n^ —0^ to vary more slowly than 
phase differences 9^^^ — 9^. (For the case of the annular 
billiard, a semiclassical argument is given in Ref. p2t.) 

Consequently, eigenphase splittings SOn'^^'^^'^ will show 
fluctuations on two scales: There will be a rapid sequence 
of peaks due to avoided crossings of regular eigenphases 
with chaotic ones and a slow modulation due to the rel- 
ative motion of regular and beach eigenphases. 

Let us conclude by summarizing those predictions that 
genuinely depend on the explicit inclusion of the beach 
layers into the five-block matrix model: 

(I) Eigenphase splitting: Contributing paths typically 
pass through all blocks. As a function of an external 
parameter, the splitting varies on two scales: a slow one 
attributed to the change of the d^^, and a rapid one 
attributed to the change of the d~\^. Consequently, one 
sees resonances of different line shapes. 

(II) Eigenphase shift: Paths contributing to the shift 
typically visit only the beach layer. The shift is much 
larger than the splitting, and it varies with df^\ on the 
slow scale only. 

None of these statements would hold for the three- 
block model, and therefore (I) and (11) can serve as a 
test of our five-block model. 
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D. Numerical Results 

We now give numerical examples of the formulas just 
presented. In particular, we will give the most quantita- 
tive and direct proof of the chaos- assisted tunneling pic- 
ture yet. Also, we will test the predictions derived from 
the explicit treatment of the beach layer in our five-block 
matrix model. Again, we consider the annular billiard at 
parameter values k = 100, a = 0.4 and i5 = 0.2. 

First of all, we have to decide on where to set the 
borders between the different blocks. As was already 
mentioned, the outer borders between the regular blocks 
and the edge blocks do not pose any problem as in both 
blocks, matrix elements are approximated by the corre- 
sponding matrix elements of the original ^-matrix. Due 
to the tunneling ridges in the region £ > k{a + S), paths 
starting from high n will most likely tunnel into this re- 
gion first. In order to include these paths, we extend the 
beach region into the regular block whenever necessary. 

The border between the beach region and the chaotic 
block is more difficult to determine. In Section [VE we 



give numerical evidence that some of the beach states' 
structure arises from trapping of classical motion near 
KAM-like regular island extending into the chaotic sea 
down to impact parameters \L\ « 0.55. This would sug- 
gest the choice of \£\ = 55 for the borders between the 
edges and the chaotic block. However, chaotic states that 
can carry transport between positive and negative an- 
gular momenta have sizeable overlap only with angular 
momentum components between £ = — 50 and £ = 50. 
Therefore, we will take the chaotic block to extend over 



angular momenta \£\ < £coe = 50. The uncertainty 
inherent in the choice of £coe affects the final results via 
the effective couplings 5^,^, Sn,'y which can therefore only 
be determined up to an overall constant of order one. 

Let us briefly discuss the magnitudes of the tunnel- 
ing amplitudes involved in the different contributions to 
the splitting. As expected, the direct n to —n tunnel- 
ing matrix element is of negligible size. At the present 
parameter values S0^g^^ ~ 1570, _7o| ^ 10^^° which is 
by many orders of magnitude smaller than the observed 
splitting SOro ^ 10" The difference between con- 
tributions S0n^"^ and 69n'^'^'^^^ is less drastic. At the 
parameter values considered, effective coupling elements 
Sn/Si^y/dn^y are usually by at least an order of magni- 
tude larger than the corresponding Sn,-y 

The dominance of recer contributions is particularly 
clear when studying the behavior of shift and splitting as 
a function of an external parameter. We show in Fig. |^ 
the shift and the splitting of the doublet n = 70 as ob- 
tained from numerical diagonalization of the ^-matrix as 
a function of the outer circle's radius R = 0.985-1.025. 
R varies over a sufficiently small interval as to leave the 
classical billiard dynamics essentially unchanged. The 
choice of R as external parameter has the advantage that 
the tunneling magnitudes given by the inner scattering 
matrix \Sn,e\ = \S^l\ remain constant, as variation of R 
affects only the outer scattering matrix S . 

Fig. |^(a) displays the real part of the shift A6n^\ 
The shift varies slowly as a function of R and is over 
long ranges well reproduced by just a single term of 



Eq. (32c). (The dashed line shows l'S'70 59/^70,59! with 




FIG. 7. (a) Eigenphase shift A6„ and (b) splitting S6n of the doublet n = 70 for different values of the radius R of the 
outer circle: exact results as obtained by numerical diagonalization (circles) and contribution predicted by the block matrix 
model that is dominant between R ~ 0.99 and R ~ 1.015 (dashed line, see text). Inset: power law dependence of splitting near 
Rp = 0.9946 (full circles) and Rp — 1.0043 (empty circles). Throughout, full lines are to guide the eye. 
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(^70,59 ~ ^70^ — ^59 extracted from S, slightly shifted to 
account for A^sg.) Over relatively large ranges of R, one 
single £ is clearly dominant. Transitions between domains 
of different dominant £ are marked by very small shifts 
due to cancelations between contributing paths. Fig. 0(b) 
shows the splitting of the same doublet. Note that the 
logarithmic scale in plot (b) ranges over twice the number 
of orders of magnitude than in (a). As predicted by the 
five-block model, the splitting is much smaller than the 
shift and shows variations on two scales. The overall, 
slow modulation is determined by the beach resonance 
DC d^Q^g and closely follows the behavior of the shift. On 
top of this modulation lies a rapid sequence of spikes 
that we attribute to quasi-crossings with eigenphases of 
the internal block. 

For the £ = 59 contribution to the splitting (dashed 
line), we used Eq. (|5^) of the next Section to estimate 
the median taken over the properties of the chaotic block 
and divided out a factor ^ 15 to make the dashed line co- 
alesce with the splitting away from the (i„^-y-resonances. 
Note that the change of dominant edge index £ is signaled 
by a strong cancelation of tunneling paths. In the inset, 
we compare the line shapes of the two types of resonances 
by plotting |^07o| as a function of | i? — i?p | in a double 
logarithmic plot near the "beach" peak {Rp = 0.9946, full 
circles) and a "chaotic" peak {Rp — 1.0043, empty cir- 
cles). Power laws with exponents —2 and —1 are obeyed 
to good precision, thus confirming the prediction of the 
five-block model. 

It is evident that the predictions of the five-block model 
serve very well to explain the data. We stress again that 
the effects just described — different line shapes and fluc- 
tuations on two parameter scales — genuinely depend on 
the role of the beach layer in tunneling processes. They 
serve as clear fingerprints of the quantum implications 
of the presence of an beach layer between phase space 
regions of classically regular and chaotic motion. 

Let us however mention that the correspondence be- 
tween shift and splitting can be less clear. For large n 
the shift can show additional modulations that do not 
appear in the splitting whenever there is a degeneracy 
with an beach mode with large £ |^l|. In the shift, this 
resonance is weighted with l^n^^ p, which favors £ near the 
tunneling ridge. In the splitting, the resonance's contri- 
bution has the weight which can become very 
small if £ is too large. This does however not contradict 
the predictions (I) and (II), it merely means that shift 
and splitting arise by coupling to different beach modes. 

For later purposes it is important to note that when- 
ever the same beach state is dominant in both shift 
and splitting, the correspondence between the shift and 
the slow modulations of the splitting can be used to 
"unfold" the splitting data from beach properties. By 
Eqs. (32c. 33a), the ratio 
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Sn,£ 



15; 



i,c\ 



r1 



(38) 



then contains only properties of the center block and can 
therefore be used to extract its "bare" quantities. 



E. Evolution of Tunneling Flux on the PSOS 

Let us recall that the scattering matrix S is the quan- 
tum analogue of the classical Poincare mapping; it consti- 
tutes a time domain-like propagator in the representation 
fixed by Eq. (^) . This makes it possible to study the evo- 
lution of " wave packets" — vectors ctQ corresponding to 
initial conditions localized in phase space — under the 
action of S. In particular, it is here of interest to follow 
the evolution of a tunneling process in phase space. 

The comparison of quantum dynamics and classical 
phase space can, in the context of the scattering ap- 
proach to quantization, conveniently be done by use of 
Wigner- and Husimi-like functions of quantum operators 
|^ , P2| . Let A be some operator that, for definiteness, 
we represent in angular momentum basis. As explained 
in detail in Ref. |22f| , A can be transformed to a function 
p^[A]{j, L) on the Poincare cell C-p by first performing 
a Wigner-transform on A and then smoothing with a 
minimal-uncertainty wave packet. One obtains 



exp ■ 



-A72 



£^£' 



+ (^ - £^ f 



- ij(£ - £') 



where {j,L) are the coordinates in C-p, k is the wave 
number, and Aj'^ is a parameter determining the shape 
of the smoothing wave packet. We choose Aj'^ = 4/k. 
In case that A ^ a ■ ot^ is the projector of a (normal- 
ized) vector a, its transform p^{j, L) constitutes a posi- 
tive semi-definite, normalized density distribution on the 
Poincare cell, the Husimi Poincare Density (HPD). By 
use of HPDs, it becomes possible to follow the phase 
space evolution of a tunneling process from one regular 
torus to its counterpart. The "dynamics" (in iteration 
count A'^ as the time variable) of such a process is visu- 
alized by projecting iterates ojn = ujq of some initial 
vector Wo onto Cp. 

Returning to the annular billiard, we consider a start- 
ing vector ujq peaked at high angular momentum n and 
calculate the Husimi densities p^[(a;jv • '^\)] of averaged 
autocorrelations 
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FIG. 8. Spread of tunneling flux over the chaotic sea; build-up of quantum amplitude in the beaches. Grayscale plots 
(arbitrary scale) of averaged autocorrelations of the (a) 50-th, (b) 100-th, and (c) 500-th iterate of a vector peaked at n = 66. 
The initial component n has been left out. Full lines indicate classical tori. 



in which the mean over 50 iterations has been performed 
in order to average out the internal dynamics of the cen- 
ter block. The tilde indicates that the n-th component 
of u) has been set to zero. (This truncation is necessary 
because the smoothing tails of the large n-th component 
would obscure all features in the nearby beach regions.) 

Fig. H depicts the flow of tunneling probability on C-p 
by showing p^[(a;Ar • a;]^)] for n = 66 and TV = 50, 100 
and 500. At these parameter values, the tunneling pe- 
riod between WG tori is 2-k/59qq ~ 10''. As predicted 
by the five-block matrix model, probability is fed from 
the starting angular momentum n into the nearby beach 
region until it reaches a value ^ \Sn,t\'^ = 5^) ^^"^ 
spreads over the chaotic sea up to a value ^ \Sn,iSifi\'^ , 
see Fig. ^(a). Oscillations between the beaches set in 
with period 27r/(5058 ~ 1000, see Fig. §(b,c). On a much 
larger time scale, probability amplitude starts to build 
up at at —n (not shown here). One clearly sees that the 
shape of the HPD is structured by the underlying classi- 
cal dynamics: in the beaches, most probability builds up 
around the chains of small KAM-like islands, whereas in 
the chaotic sea, the center island and its satellites are not 
penetrated. Also, the regions around the satellite islands 
and the homoclinic tangles between them are filled only 
weakly. 

Chaotic phase space can also be filled in a different 
manner, depending on the phases On, de and 9^ involved 
in the tunneling process. In Fig. ^ we present the case 
of a close degeneracy between 9n and one of the 9^. We 
show (a; AT • uj]^) for N = 4000 and the starting vector ujq 
peaked at n — 65. In this case, there is high probabil- 
ity amplitude in the sticking regions around the center 
satellite islands. 



Finally, we present in Fig. |10| HPDs of nine eigenvec- 
tors of 5* at fc = 100, a = 0.4 and S = 0.2 (k = 100, 
one should note, is not an eigenenergy of the annu- 
lar billiard). By the Weyl formula these vectors 
would, when quantized at a close-by energy, correspond 
to the ^ 2000-th exited states. We show grayscale plots 
of HPDs with steps in grayscale corresponding to 
equidistant probability contour lines. The overall scales 
vary with each sub-plot. These HPDs at hand, the spread 
of tunneling amplitude can now be understood in terms 
of participating eigenvectors. For example, the doublet 
depicted in (b) and (c) is the beach doublet involved in 
the tunneling process of Fig. ^ Likewise, the vector (f) 
peaked around the center satellite islands is the nearly 
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FIG. 9. Spread of tunneling flux over the chaotic sea: case 
of a direct degeneracy between the regular doublet and an 
internal state. Grayscale plot (arbitrary scale) of averaged 
autocorrelation of the 4000-th iterate of a vector peaked at 
n = 65. 
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FIG. 10. Poincare Husimi Distributions of selected eigenvectors of the annular billiard at fc = 100, a = 0.4 and 5 = 0.2. (a) 
Regular high-angular momentum vector, (b,c) doublet of beach vectors, (d,e) "chaotic" vectors, (f ) vector in the sticking region 
around the satellite islands, (g) regular vector on the period-6 satellite islands, (h) regular vector residing on the main island, 
and (i) vector scarred by the unstable period- 1 fixed point and its homoclinic crossings (dots depict the stable and unstable 
manifolds) . 



degenerate one in the tunneling process shown in Fig. 
Evidently, their shape forms the spread of tunneling 
probabihty on the Poincare cell. 

We note that similar figures for much higher wavenum- 
ber k — 600, corresponding to the ~ 75000-th excited 
states, can be found in Rcf. 12^. 



F. How Important is Chaoticity? 

Let us now discuss a numerical study similar to that 
performed in the original work by Bohigas et al. in 
which regular level splittings are calculated as a function 
of eccentricity 5 at constant a + 5. It is important to 
note that increasing 5 has two effects: classical motion 
in the inner layer \L\ < a + 5 becomes chaotic, and si- 
multaneously tunneling rates from the regular torus to 
the chaotic layer are enhanced. A priori^ it is not clear 
which one of these effects governs the rate at which the 
splittings change, but a quantitative estimate of either 
effect has now become possible. 



We calculated the splittings of high-angular momen- 
tum modes for 5 — 0.03-0.2 and a + 5 — 0.6 at A; = 60. 
Fig. |ll| displays the eigenphase splitting \56n\ for n = 39 
(n/fc = 0.65). Exact splittings (full lines) increase over 
eight orders of magnitude and roughly follow an expo- 
nential increase with 5. A description of the data in 
terms of our block-matrix model that reproduces all the 
fine details might be a difficult task — even with ex- 
act S'-matrix elements at hand — as the model relies on 
classical information to select the borders between the 
different blocks and assumes that, apart from the beach 
layers, no significant phase space structure is present. 
Therefore, the block-matrix model would have to be ad- 
justed to the varying classical dynamics as 5 changes, and 
the effect of remaining structure at lower 5 might have to 
be taken into account with the introduction of different 
blocks. However, the now-familiar slow modulations in 
Fig. |ll| point to the effect of beach layer states mediating 
the tunneling flux — however complicated the internal 
structure might be. Indeed, we find that over the range 
5 — 0.07-0.16 the tunneling processes are mediated by 
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two beach states peaked around t = 29, with ijk = 0.48 
weh inside the non-integrable regime. In Fig. |ll|, we have 
plotted |S'„/p as a dotted Hne (with arbitrary offset) to 
give a rough estimate of the change of regular-to-beach 
tunneling elements via these particular beach states. 




FIG. 11. Eigenphase splitting \ &Qn\ of doublet n = 39 as 
a function of eccentricity 5: exact splittings (full line), in- 
crease of torus-to-beach tunneling matrix elements [i = 29 
dotted line, with arbitrary overall factor), and effect of reg- 
ular-to-beach phase denominators (dashed line, with arbi- 
trary overall factor). For 5 > 0.15: tunneling via £ = 35 
(dashed-dotted line). 

Taking into account also the effect of beach denomi- 
nators, we have plotted \Sn,ildn,t\^ (dashed line). This 
estimate is multiplied with an overall factor 2 • 10~^ to 
make the line coincide with the exact data at small b. 
The dashed line already gives a fair reproduction of the 
data. It misses only the sharp peaks due to resonances 
with the chaotic states, the depression of (50„ between the 
d~ £ peaks due to destructive interference, and the change 
of coupling strength of the beach state to the chaotic 
center states. We conclude that between (5 — 0.07 and 
(5 = 0.2, the splitting is predominantly determined by the 
change of beach parameters, and that the change of in- 
ternal coupling between beach and chaotic sea accounts 
only for a factor of the order 10 (difference between the 
dashed line and the exact data a,t 6 = 0.15). For larger 
6, different beach doublets take over, but the basic struc- 
ture is preserved. The dashed-dotted line displays the 
recer-contributions for £ = 35 and £coe = 20. 

At ^-values below 0.07, numerical precision does not 
allow us to calculate the splitting directly. We can how- 
ever get an impression by looking at the splitting of states 
supported by quasi-integrable structures at small S. Pre- 
sumably, these states will mediate the tunneling of high- 



angular momentum doublets. In Fig. |lj, we show the 
splitting of the doublet predominantly peaked at £ = 27 
{£/k — 0.45). The splitting shows resonance peaks be- 
low 6 w 0.075 and then flattens out. This behavior can 
be understood by looking at the HPDs of the states in- 
volved. In Fig. |l3| we depict one partner of the tunneling 
doublet and its resonant state at (a-c) S = 0.0402, (d~f) 
S = 0.0687, and (g-h) S — 0.09. In the corresponding 
classical Poincare cells, we have started trajectories from 
initial conditions (tt, L) with L > only to indicate the 
classical inhibition of transport. We can make two inter- 
esting observations: First of all, the resonant tunneling 
process is mediated by regular states residing on the inner 
island. Resonant tunneling via the center island remains 
the dominant mechanism, even when classical transport 
from positive to negative L becomes allowed (d-f). Here, 
tunneling between a chaotic doublet is mediated by a reg- 
ular state. Secondly, the outer doublet supported by the 
KAM-like tori at small 6 evolves into a doublet of states 
scarred near the unstable periodic orbit and stretched 
along its stable and unstable manifolds. By quantum 
localization effects, the doublet structure persists — de- 
spite of the seemingly chaotic classical motion (see also 
p2t). Tunneling between the scarred doublet is direct, 
as indications of resonances are absent in the splitting 
beyond 5 = 0.075. 
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FIG. 12. Eigenphase splitting of a doublet peaked around 
£ = 27 at fc = 60 as a function of 5. Here, resonances arise 
from avoided crossings with states residing on the center is- 
land. 

We are led to the conclusion that the enhancement of 
tunneling rates between symmetry-related phase space 
objects A and TA by resonance with quantum states 
supported by an intervening phase-space structure B is 
only very loosely related to the chaoticity of B, but rather 
depends on the topological character of B. In B, it must 
merely be possible to traverse phase space distance in 
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FIG. 13. Classical phase space and Husimi Poincare distributions of the tunneling doublet of Fig. |l^ and its resonant states 
at (a-c) S = 0.0402, (d-f) 6 = 0.0687 and (g-h) S = 0.09. The doublet resides on KAM-Uke tori for small S and dissolves in 
the chaotic sea for large S. Note that the tunneling at 5 = 0.0687 takes place between a chaotic doublet and is mediated by a 
regular state. At S = 0.09, the dominant tunneling process is direct. 

classically allowed steps. It might therefore be more 
appropriate to refer to "transport-assisted tunneling," 
a phenomenon of a more general class than the chaos- 
assisted tunneling one. 



Starting from Eq. (32c) we introduce a number of no- 
tational simplifications. For a given n, we write phases 
with respect to 6'i°'' (that is, set = 0), define 

= 2sin(0^/2) and collect the coupling of ±n to a 
chaotic state 7 via the edge into an effective overlap 



V. STATISTICAL ANALYSIS OF EIGENPHASE 
SPLITTINGS 

A. Asymptotic Behavior of the Splitting 
Distribution 

This Section is devoted to the distribution of level split- 
tings. We determine the asymptotic large-(50„ behavior 
of the splitting distribution and find "typical" splitting 
values by calculating the median of the distribution ob- 
tained by extrapolation of the asymptotic behavior to- 
wards smaller splittings. We assume that recer contribu- 
tions dominate and that properties of the beach blocks 
vary slowly. It is straightforward to apply the calculation 
to the case of rcr contributions as well. 
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which leads to 



1 (recer) 
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} , (39) 



Nj denotes the dimension of the center block. In the 
sequel we drop the subscript n. 

To devise a statistical treatment for the center block, 
we make the following assumptions concerning the dis- 
tribution of the chaotic eigenphases 9 — {9^}, and the 
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effective overlaps v = {v^}, 7 = 1,..., Nj. We assume 
that 

(1) there is no correlation between the overlaps and 
the eigenphases, 

(2) eigenphases O-j are real, ranging from [— 7r,7r], and 
the joint distribution function P{6) of the eigenphases 
either (a) is Poissonian, i.e. the eigenphases are uncorre- 
lated, or (b) has the property that degeneracies of eigen- 
phases 0y are suppressed (as is the case in Dyson random 
matrix ensembles), 

(3) the overlaps are mutually independent Gaussian 
random variables with zero mean and variance cr^ . a will 
be fixed in terms of ^-matrix elements in the following. 

The joint probability distribution of the 6j and v-y is 
then given by 



1 

p{e,v) = Pie)l[ — 



7=1 



We can write the probability density of 66 in the form 



P{S9) 
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(40) 



where P{66\ 6) is the conditional probability of 66 given 
6, and the integral is performed over [— tt, Tr]''^^. For fixed 
6, the v-y/x-y are mutually independent Gaussian random 
variables with variances (cr/7)^, and thus P{59\9) is a 
Gaussian of variance G^ri{6), where rj{6) is given by 
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7 — 1 ' 



To perform the integration over 9 in Eq. (^0|), we in- 
troduce 77(0) as an additional integration variable by re- 
writing P{S9) as 
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with 
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(41) 
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Note that, by virtue of the Central Limit Theorem, al- 
most all overlap distributions will give rise to a Gaussian 
form of P(b6\ 6) in Eq. (|4^). 

We continue by examining the large-77 asymptotic be- 
havior of P{ri) for the two cases listed in condition (2). 

(a) If the 9j obey a Poissonian distribution, then 77 is 
a sum over independent identically distributed variables 
r]^ , each of which is distributed as 



P(,,).^£d.. 







1/4] 



for 
for 



»77 



> 



'77 < 



4 ' 
1 

4 • 



Note that for large % this asymptotically behaves like 
P{rij) ^ l/27rry7^^. In order to obtain the distribution 
P{t]) of the sum, we evaluate the characteristic function 
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(see H, Eqs. 3.383.4 and 9.236.1]) where erfc(a;) denotes 
the complementary error function. We use the fact that 
the characteristic function of a sum of Ny independent 
random variables is the product of all their characteristic 
functions, to get 



duj erfc 



The structure of the 77 ^ 00 tail is determined by the 
non-analytic behavior of Pi{i-u) at the origin. Thus we 
can approximate P{i]) in that regime by expanding in a 
power series in lu, 
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The leading (non-analytic) square root term is propor- 
tional to Nj. Thus the frequency dependence scales with 
N^, and the resulting distribution P{ri) scales asymptot- 
ically as P{rjlN'^). We thereby obtain the asymp- 
totic distribution of 77, 
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(b) Next we will evaluate P{ri) for a general eigen- 
phase distribution in which the occurrence of eigenphase 
degeneracies is suppressed Jk5[. In order to re- write the 
distribution function Eq. (|42| ) in terms of products over 
the eigenvalues, rather than sums over them, we intro- 
duce the integration parameter a = by writing 

P{v) = f_Aajdep{e)5(^-Y,^ 5(^-\{x^ . 

In the first (5-function, we can now substitute 

E ^ ^ ^ E II = ^ E '^7 ' 

^7 7 y'^y 7 

where we have introduced the notation X-y — Jl7'5^7 ^7'- 
By changing the integration variable a ^ oij ^Jt] we 
can extract the explicit r/-dependence from the first 5- 
function and arrive at 



P(r7) 



o3/2 



dep[e) 



da^ 1 



nl/2 



7 / 



n 



(44) 
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By inspection of the first (^-function, and recalling that 
< 1, we see that contributions to the integral can 

1/2 

only arise from the range \a\ < N^' . However, for finite 
Of, the second (5-function is in the limit rj oo given by 




J2 H^i) , 



(45) 



provided all the distinct. Hence, the second S- 

function is asymptotically independent of a, and the a- 
integration can be performed explicitly, using 



1/2 



(46) 



Finally, we substitute the form (^5|) of the second 5- 
function and the result (^6|) of the a-integration into 
Eq. (H). Noting that if 0, then Xy/ = for aU 

7' 7^ 7, and that S{9'y) = ^(x^), we arrive at the distri- 
bution 



27r?73/2 



for 



V 



(47) 



It is remarkable that the distributions P{r]) in Eqs. ( p^ ) 
and (|4^) display the same asymptotic power-law depen- 
dence. 

Finally, we note that, by virtue of the Gaussian form of 
the integrand in Eq. ( ^ ) , the integral depends primarily 
on large rj > 69^ /a^. To extract the asymptotic behavior 
of P{59), we only need the asymptotic form of P{ri) for 
large 77. Inserting (^^ into Eq. ( ^ ) gives 



P{S9) - 



dri 



-5e^/2j;CT^ 



(27r)3/2^2^ 7rV2^(56l2 



(48) 



for large 69. This is the large-69 tail of a Cauchy distri- 
bution, in accordance to the prediction of Leyvraz and 
Ullmo d]. We stress that we have derived this result for 
the case of a Poissonian eigenphase distribution and, in a 
second derivation, without assuming any explicit form of 
the joint eigenphase distribution function. In the latter 
case, we only had to make the assumption that the joint 
distribution function vanishes whenever two eigenphases 
approach each other. Our derivation is therefore more 
general than the one given in |^]. 

In a broader context, it is also interesting to mention 
that the calculation (b) can be generalized to the case of 
a distribution 



PAS9) = Jdep{e) sL-Y, 



7 ^7 



(49) 



which corresponds to the distribution resulting from a 
sum over paths that contains a i^-th power of the phase 
denominator. The procedure is analogous to the case 



v = 1 just described, but for a substitution a 

and one arrives at P^{r]) = (2Trv)-^ N^/ 7]-'^^"+^^^^" . 

This leads to a la.ige-69 splitting distribution 



P{S9) (X 



1 



59 (i^+i)/!' 



(50) 



Let us return to Eq. (p8|). Having integrated out the 
eigenphase dependence of P{S9), we are left with the de- 
termination of the variance of the effective overlaps 



_2 ^ 



It is given by 



2sin(0,/2)2sin(0,V2) 



We assume the phases of the Si,^ to be arbitrary and 
uncorrelated, assume the 9i to be real, and absorb all 
phases into a random phase factor exp(i(/)£^f/,7). Using 



we can write 



E 



E 



\Sn/Sn,i'Sl^^Sl'^^\ COS (j>i^i'^^ 

1 2sin(0,/2) sin[(^?,, -0i"^)/2] 
8 sin^ (6*^/2) sm'^{9i-/2) ' 



(51) 



using the fact that the (peji^y are uncorrelated and equi- 
distributed on the interval [— tt, tt]. Let us now write 
|5,,7|2 = N-^Si,cHl^, ^l^e^e l^^.cP = E7 1^^.71' is 
the total coupling of the £ state to the chaotic block, 
see also Eq. (p8|), and where the ^^^^ are indepen- 
dent Gaussian variables with unit variance. Using that 
( 7^1' 7 ) = (2 + Se,e') we find a statistical enhancement 
of the diagonal {£ = I') terms. More importantly, the 
sum is dominated by the terms with the smallest phase 
denominators. Consequently, we can neglect the non- 
diagonal terms and write 



8^2 



E 



sm{9il2) 



(52) 



B. Median Splittings 

As we have just shown, the splitting distribution be- 
haves asymptotically like P{69) ~ 69 and it is well 
known that the mean of a Cauchy distribution does not 
exist. Therefore, a "typical" value for the level splittings 
must be obtained otherwise. We propose to consider the 
median \d9\M of \S9\ defined by 

diS9) pm ^ \ . 
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The factor two on the left hand side arises from the fact 
that we integrate over positive 69 only. By extrapolation 
of the asymptotic form of P{S9) as given by Eq. (^) 
towards smaller \S9\, we find 



M 



Inserting the variance of Eq. ( 
finally get for the median splitting 



(53) 

l5|) into Eq. (|5|) we 



M,n ~ ~ 




.C 



sin[(0,-^^i°')/2] 



1/2 



(54) 



where we have inserted the index for the n-dependence 



again, as well as the phase 9n^\ Formula (54) for the me- 
dian splittings estimates the enhancements of tunneling 
splittings due to chaos-assisted processes and constitutes 
one of the central results of this work. Note that all 
quantities appearing in Eq. (|^) are defined in terms if 
the original ^-matrix, and the most a direct and quan- 
titative check of the chaos- assisted tunneling picture yet 
becomes possible. 

We note that Eq. ( ^4| ) can be used to recast Eq. ^ ) 
for the splitting distribution in a numerically more con- 
venient form 



\69\ 



M.n 



A 591 



(55) 



We now turn to the discussion of the approxima- 
tions made in the derivation of the central results 
Eqs. (^,^^5|) for the splitting distributions and the 
median splittings. There are four sources of error, (i) 
The estimate ( p^ , p9[ ) for the variance is correct only 
within an order of magnitude due to the ambiguity of 
^COE (see discussion after Eq. (|29|)). Up to now, an a 
priori determination of ^coe is not possible. Note, how- 
ever, that the size of the center block does not enter 
in the expressions, (ii) The effect of imaginary parts of 
the 9~f is not included in our calculation. By unitarity 
of the block-transformed matrix, 3m{0^} ~ |5£_-yp, 
which in the splitting distribution introduces a cutoff of 
the Cauchy-like tail at \59\ ~ \Sn,i/ s\n{9i/2)\'^ . The 
resulting relative correction of the median splitting is 
0{\Si,c\'^). (iii) Extrapolation of the asymptotic tail to- 
wards smaller 59 is another source of error of 0(1). (For 
example, the median calculated from an exact Cauchy 
distribution l/7r(l 4- x^) is 1 + \/2 w 2.4, whereas the 
median estimated by integrating over its tail I/ttx^ is 
4/7r K, 1.3.) (iv) Our five-block model neglects the ef- 
fect of transport barriers other than the one separating 
the beach from the center of the chaotic block. Further 
transport barriers lead to the inhibition of tunneling flux 
and thereby decrease the splitting. 



We conclude that (54) reproduces the exact median 
splittings only up to a factor of the order one. If the 



neglect of remnant phase space structure is the domi- 
nant source of error, then Eq. (^^ gives an over-estimate. 
However, the error is expected to be independent of n, 
and we can correct for it by introducing an overall factor 
c that we extract from the numerical data. The formula 
Eq. ( p5| ) for the splitting distribution function has to be 
corrected correspondingly. 

We finally return to the issue of the two different 
representations Eq. (p3|) and Eq. (24) for the split- 



ting that differ by taking either imaginary parts of 



e-x^{—iN9n'')[S'^]n,-n or absolute values of [S 
a statistical treatment, these two approaches give slightly 
different results, because in the average over the random 
phases ip^/'j, one obtains (le^^^ ^ T^I) = 1 after taking ab- 
solute values, as opposed to {cos'^{(j)i^e,f)) — 1/2 after 
taking imaginary parts. The median splittings derived 
from Eq. (|2^) would therefore be twice the splittings pre- 
dicted in Eq. ( 54 ) . This explains why a formula given by 
us earlier [ p^ Eq. (8)] differs by a factor two from the 
one in Eq. (d4). 



ATI 



In 



C. Numerical Results 

This Section is concluded by a presentation of numer- 
ical data for the eigenphase splitting and its distribution 
for the annular billiard. We choose parameter values 
k — 100, a — 0.4, 5 = 0.2 and vary the outer radius 
R over 490 values between R = 0.985 and R = 0.1035. 
Recall that changing R leaves the \Sn,e\ constant and 
changes only the eigenphase configuration and the cou- 
plings Se.j. 



O 




FIG. 14. Median splittings \S6\M,n as a function of angu- 
lar momentum (logarithmic plot). Exact splittings are ob- 
tained from numerical diagonalization (full circles), calcula- 
tion of |[S^]n,~n| for large A*' (empty circles), and estimates 
are taken from the median formula ( p^ and corrected by 
c = 1/6 (dashed line). 
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Fig. |4| depicts the median values of eigenphase split- 
tings of doublets peaked at angular momenta n = 
65, . . . , 80 as a function of n. Full circles represent me- 
dian splittings |(50|A/,n as obtained from numerical di- 
agonalization. Since the diagonalization routine cannot 
differentiate between eigenvalues which are closer than 
~ 10~^^, splittings beyond n — 75 could not be resolved 



directly. Instead, they were extracted by use of Eq. (24), 
that is, by numerical calculation of 2[S'''^]„__„/A'' for large 
N ^ \Sdn\~^- Empty circles represent median splittings 
thus obtained from 30 configurations with R ranging from 
i?=ltoi?=1.3. We have taken the edge region to 
extend over angular momenta £ = 56, . . . , 64 and have 
chosen -^coe = 50. To account for the over-estimate of 
the splitting by Eq. (|5j), theoretical predictions are mul- 
tiplied by an overall factor c « 1/6. The dashed line 
shows the resulting approximation for the median split- 
tings. Apart from the factor c, the formula (Q) is in 
good agreement with the exact median splittings. 



tion of width c. For the figure, 750 exact splittings were 
transformed to reduced ones (see inset) and collected in a 
histogram with log-binning (main figure, solid line). This 
is compared to the Cauchy distribution 2cn^^ /{c^ + SB"^) 
(main figure, dashed line). The agreement is very good. 

Finally, we show in Fig. |l^ the distribution function 
P{\S9n\/\S9\M,ri) of exact splittings divided by the me- 
dian splittings displayed in Fig. ^ The actual splittings 
display a power law P(|(50|) ^ |(5^^|^'^/^ (dashed line). 
This can be understood by realizing that the variation 
in R is sufficiently large to average not only over avoided 
crossings between chaotic and regular eigenphases, but 
also over avoided crossings between regular and beach 
eigenphases. These avoided crossings, however, appear 



in the sum (32c) with a squared phase denominator d^^, 
and we have argued in Eq. ( |50| ) that the distribution 
generated by such contributions displays a power-law de- 
cay with exponent —3/2. The exponent —3/2 conforms 
to the findings of Leyvraz and UUmo who studied 
chaos- assisted tunneling in the presence of a imperfect 
transport barrier in the chaotic sea. 




FIG. 15. Distribution function P{\56\) for the "reduced" 
splitting \S6\ (solid line), compared to the predicted Cauchy 
distribution (dashed line), double logarithmic plot. Inset: "re- 
duced" splittings obtained after unfolding the modulations 
due to the beach layers. 

Let us turn to the splitting statistics. If one is inter- 
ested in the fluctuations due to changes in the chaotic dy- 
namics, one first has to discard the slow modulation due 
to the change of beach layer properties. We do so by con- 
sidering the "reduced" splitting 69 n of Eq. (38) for those 
values of R at which one single £ is dominant in both shift 
and splitting (for n = 65, . . . , 67, and here £ = n — 7). 
We find that the median of 59n is approximately equal to 
c, independent of n (not shown). Fig. 15 confirms that 




the distribution P(|(56'|) falls off like a Cauchy distribu- 



-4-2 2 



FIG. 16. Distribution of original eigenphase splittings 
\59n\, for each n divided by the median \S6\M,n (sohd line), 
double logarithmic plot. Dashed line: comparison to a S9~^^^ 
power-law decay, prefactor fitted to the data. 



VI. DISCUSSION 

A. Possible Experimental Realizations 

Even though the occurrence of chaos-assisted tunneling 
should be a very general phenomenon, an experimental 
proof of the effect has not yet been given. The main diffi- 
culty might not be to measure the effect, but to recognize 
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it. As long as little is known about tunneling in multi- 
dimensional mixed systems in general, it will be difficult 
to separate out the different contributions to the tunnel- 
ing rates and to identify the effects of classical transport. 
It is our strong suspicion that, as soon as qualitative 
theories for experimental systems are developed, chaos- 
assisted tunneling will turn out to be a frequent effect in 
the splitting of dynamical tunneling doublets. 



1. Superconducting Microwave Cavities 

It has been argued in this work that the annular bil- 
liard serves as an excellent paradigm for chaos-assisted 
tunneling. An experimental realization of it is presently 
investigated by the Darmstadt group of Richter et al. . 
Performing resonance measurements on a superconduct- 
ing niobium microwave cavity, the Darmstadt group has 
extracted high-quality spectra in the frequency range 0- 
20 GHz, corresponding to fc = 0-50 in our units, but 
experimental accuracy does not yet allow a resolution of 
the splittings of high-angular momentum doublets. How- 
ever, it might be just as interesting to measure the energy 
splittings of beach doublets, as these splittings are also 
chaos- assisted, but by orders of magnitudes larger than 
those of the regular doublets. 



2. Atomic Systems 

Atomic systems have served as paradigms of many pre- 
dictions in quantum chaos, and there are some atoms in 
which an observation of chaos-assisted tunneling might 
be conceivable. Hydrogen in a weak magnetic field ||47| ] 
may be such a case. In the weak field limit, the classi- 
cal system has symmetry-connected regular islands cor- 
responding to low-angular momentum motion along the 
field axis on either side of the hydrogen core. These is- 
lands are separated by a chaotic sea, and dynamical tun- 
neling between the corresponding low-angular momen- 
tum quantum states should therefore be enhanced by 
chaos- assisted processes. It is however not clear, whether 
these splittings are large enough to be experimentally 
accessible. Another much-studied system in the field of 
quantum chaos, the quantum kicked rotor, has recently 
been experimentally realized by Moore et al. using 
ultracold sodium atoms in pulsed, near-resonant light. 
Rotor systems have been considered by several groups in 
studies of dynamical tunneling in the presence of chaos 
pl| , ^ , and a link between theory and experiment might 
soon become possible. Again, the experiment might still 
be far from the required degree of accuracy. 



3. Open Systems 

Finally, chaos-assisted processes can enhance not only 
tunneling oscillations, but also the decay of regular modes 
in a mixed, open system in which the dominant cou- 
pling to the continuum is mediated by states residing on 
the chaotic layer. Experimental realizations of such sys- 
tems were studied by Nockel et al. and others |Q who 
considered the Q-spoiling of whispering gallery modes 
in deformed lasing droplets. In a recent work, Hacken- 
broich and Nockel [ pl| also considered mixed systems in 
which the direct coupling of regular modes to the contin- 
uum is suppressed, but where chaotic states have sizable 
coupling to the continuum. Regular modes may then 
decay via a multi-step process of type regular- chaotic- 
continuum. Their results were motivated by a study of 
a modified version of the annular billiard, in which the 
outer circle is replaced by a mirror, and the billiard is 
supposed to have higher optical density that the exte- 
rior region outside the mirror. It was found that chaos- 
assisted decay can lead a dramatic enhancement of level 
widths. 



B. Discussion and Conclusions 

Having discussed the experimental perspectives, one is 
immediately led to the question of the general applicabil- 
ity of the method and the block-matrix model proposed 
in this work. 

It is clear that the scattering problem must be solved 
separately for each system under consideration. In a gen- 
eral system, it might be very difficult to formulate the S- 
matrix and to find a basis in which the S'-matrix elements 
Sn,n and S-n,-n Corresponding to motion on the tunnel- 
ing tori are sufficiently close to unity. It must however be 
said in favor of the scattering approach that the difficulty 
of finding an EBK-quantization scheme in non-separable 
systems is by no means inherent to the scattering ap- 
proach, but presently poses one of the most serious prob- 

(In fact, it is 



lems of semiclassical theory in general 1 52 



one of the fortunate aspects of the annular billiard that 
the angular momentum basis is semiclassically diagonal 
in the region of regular motion.) 

In situations where S±n.±n are not sufficiently close to 
unity, Eqs. (^,^^ may still work well, provided that suf- 
ficient knowledge of eigenvector structure is available (see 
1^ for an application to the case of rough billiards |5^ ] 
and other systems). However, calculation of the splitting 
by summation over long paths from n to — n relies on a 
sufficient localization of |n) at the n-th component. 

Clearly, the block matrix model used in the summa- 
tion over paths must be adopted to the specific transport 
situation encountered. In case of structure other than 
the beach layers, additional blocks must be introduced. 
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This does however not lead to problems in the summa- 
tion Eq. as long as the outermost tunneling element 
is smaller than any of the internal coupling elements. At 
present, there is no a priori method to determine the 
border indices of neighboring blocks inside the chaotic 
layer. However, use of classical information will warrant 
correct results to within an order of magnitude. 

We note that our treatment is not limited to the case 
of an S'-matrix symmetry Sn,m = S-n,-m and could 
easily be extended to non-symmetric systems (for exam- 
ple, an annular billiard with the inner circle replaced by 
some non-symmetric shape) , or even to the case of tun- 
neling at an accidental degeneracy between two eigen- 
phases. In the summation over paths, we merely re- 
quire that the initial and the final diagonal S'-matrix el- 
ements are equal, Sni,ni = Snf,nf Note, however, that 
in the non-symmetric case the contributions {ni\S^\ni) 
and {nf\S^\nf) in Eq. (Q) will in general not cancel. 
Their difference is then likely to dominate the splitting. 

It is one of our main results to point out the impor- 
tance of the beach layer to the chaos-assisted tunneling 
phenomenon. The appearance of classically chaotic, but 
not too unstable regions around regular islands is generic 
in mixed systems. Such regions should always support 
states if the mixing with the rest of phase space is suffi- 
ciently slow (or if energy is sufhciently low). However, it 
must be checked whether some of the importance of the 
beach layer should actually be attributed to the tunnel- 
ing ridge that favors tunneling processes into the beach 
region (see Fig. |^). Such a test is given in our version of 
the Bohigas numerical experiment, see Section IVF. We 



rather insensitive to the rate of classical flux connecting 
the opposite beach regions — as long as there exists a 
classically allowed path between them. When changing i5 
in Section [VF, most of the tunneling enhancement was 



have verified that, at small 5, there are no visible tun- 
neling ridges, and the tunneling amplitudes decay mono- 
tonically away from the diagonal. Nevertheless the beach 
region still governs the behavior of the eigenphase split- 
tings (see Fig. |l^). It must however be noted that, in this 
case, the correspondence between slow splitting modula- 
tions and the shift breaks down. The shift is then deter- 
mined by paths n ^ n — 1 ^ n instead of paths leading 
to the beach and back to n. 

It is another interesting point that the statistical re- 
sults found in Section ^ are independent of the explicit 
joint distribution function of eigenphases, but can be de- 
rived under rather general assumptions. We merely re- 
quire that the joint distribution of eigenphases vanishes 
at eigenphase degeneracies. 

This generality is extended even further by the ob- 
servation that enhancement of tunneling can also ap- 
pear with help of regular states. Indeed, in Section IVF 
we have even observed the case of tunneling between a 
chaotic doublet via a resonant regular state. This should 
serve as a reminder that only the phase space topology 
determines the occurrence of tunneling, not its regular- 
ity or chaoticity, and that chaos-assisted tunneling is, in 
fact, a more general phenomenon of transport-assisted 
tunneling. Additionally, the tunneling rate seems to be 



related to the change of tunneling properties between the 
torus and the beach region. Progressively rapid classical 
propagation across the chaotic layer was related to tun- 
neling enhancement of only one order of magnitude — 
out of five orders of magnitude in total [5 = 0.07-0.15). 

Our study of chaos-assisted tunneling has led to the 
most quantitative treatment of the phenomenon to date. 
At the same time, some challenging problems have been 
encountered. For example, we have seen that tunnel- 
ing can occur between doublets localized on "soft" phase 
space structures such as the beach regions or scarring 
periodic orbits. For these states, transport from one 
phase space structure to its symmetry-related partner is 
classically allowed, but quantum mechanically forbidden. 
Apart from the intriguing question about the quantum- 
mechanical localization mechanism giving rise to these 
states, their doublet structure introduces additional com- 
plications. For example, some of the doublets tunnel via 
resonant processes, while others tunnel directly. A quan- 
titative treatment would certainly be desirable. 



C. Summary 

We studied dynamical tunneling between symmetry- 
related phase space tori that are separated by a chaotic 
region. Using scattering theory, we introduced a unitary 
matrix S that constitutes the quantum analogue of the 
classical Poincare map. By expressing eigenphase split- 
tings and shifts in terms of matrix elements of high iter- 
ates of S, we related these quantities to paths in phase 
space. While paths contributing to the splitting connect 
the two tunneling tori, paths that contribute to the shift 
lead from a tunneling torus back to itself, leaving the 
torus at least once. We performed the summation over 
paths within a block-matrix approximation, allocating 
different blocks to the two regular regions, the chaotic 
sea, and the two intervening beach layers. Within this 
approximation, we derived analytic expressions for the 
contributions to the tunneling properties. Explicit inclu- 
sion of the beach blocks enabled us to predict a number 
of new effects that could be verified for the case of the an- 
nular billiard. (I) As a function of an external parameter 
the splitting varies on two scales: a rapid one attributed 
to resonance denominators of regular and chaotic states, 
and a slow one attributed to (squared) resonance denom- 
inators between regular states and beach states. This 
diversity of scales is also observed in statistical quan- 
tities, e.g the distribution function P{56) of eigenphase 
splittings 59. When averaging over a sufficiently small 
range of a system parameter (such that beach proper- 
ties remain effectively constant), the splittings are dis- 
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tributed with a Cauchy tail P{59) ~ S9~'^. When av- 
eraging over a large parameter range (such that beach 
resonances occur), the squared resonance denominators 
lead to a P{59) ~ 69~^/'^ power-law behavior. (II) Typi- 
cally, the shift varies on the slow scale only and is much 
larger than the splitting. 

Analytical formulas at hand, we could also asses the 
relative importance of tunneling amplitudes and clas- 
sical transport properties within the chaotic sea. As 
the annular billiard's eccentricity is increased, most of 
the enhancement of tunneling rates can be attributed to 
the tunneling amplitudes and resonances between regular 
tori and the beach regions. Progressively faster classical 
transport within the chaotic sea was found to play a mi- 
nor role in the splitting enhancement. 

Finally, we derived the asymptotic form P{59) ~ 156*"^ 
of the splitting distribution's \aige-59 tail (average over 
a small parameter range). In this calculation, no ex- 
plicit assumption about the form of the joint distribu- 
tion function of chaotic eigenphases was made; it was 
merely required that the distribution is either Poissonian 
or vanishes for degenerate eigenphases. In order to give 
"typical" splitting values, we calculated the median split- 
ting by averaging over the properties of the chaotic sea. 
Apart from an over-estimate by an overall factor ~ 5, the 
predicted values for the median splittings closely follow 
the numerical results over many orders of magnitude. 
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a grant from MINERVA. 



APPENDIX A: SUMMATION OVER PATHS IN 
THE SPLITTING FORMULAS 

In this Appendix we will derive formula ( ^l]) that con- 
tains the contributions of the different families of paths 
to the full sum over paths 



^N 



N-1 

E 

{n — n} i=l 



n ^ 



Ai.Ai 



leading from n to —n. Recall that 'P,^_„ is related to the 
splitting of the doublet 59^ by Eq. @j. As N ^ K/d9n is 
taken to be large, we need only collect the leading-order 
contribution in N . 

Let us consider the general case in which the S'-matrix 
and block-diagonalized into any number of blocks S^^'^\ 
i,j = 0,...,K. For the summation over paths, we 
merely require that the outermost couplings sf'y and 

S^'y"^ be much smaller than entries of the internal cou- 
pling matrices. Suppose we want to collect the contri- 
butions from paths with M + 1 steps that start from n 



and pass through the M intermediate diagonal blocks 
S^'^^'^\S^'^^'^\...,S^'"^''''> before arriving at -n, and 
let us for the moment neglect repetitions in the block in- 
dices. It is equivalent to sum over all such paths in the 
block-tridiagonal matrix 
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Ci D2 








c 



V 



Cm- 



Dm 
Cm 



C'm 
Do J 



(Al) 



where each of the diagonal blocks Di — u^y^ 
pled to its neighbor by the coupling block C, 



S\ i Oa.a' is con- 



s' 



{Li+l) 
A, A' 



C is a short notation for the other coupling block, C[ — 



S 



{i+lA) 
A'. A 



(which need however not be the transpose of C;, 



as S is in general not symmetric). Do = S'n,n = S-n,-n 
contains the diagonal S-matrix element associated with 
the tunneling tori. 

In addition to summations over the matrix indices 
of each block, we have to sum over the staying times 
No, ... , Nm inside the diagonal blocks Do, . . . ,Dm, re- 
spectively. Explicitly, we have to perform the sum 
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{A'o,---A'm} 



Do""" Co D^' ■ ■ ■ D^l' Cm i?Jv-Jv„-JVi...-Jv,,-M 
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J2 Co{D,/Do, 



{Dm/Do)^-' Cm 



{No,-Nm} 



where the summations go over Nm — 0, . . ■ , N — M, 
Nm-1 = 0, . . . , iV - Nm - M, and so on, up to No = 
N — Nm — Nm-i ~ . . . ~ Ni ~ M. Summing over A^'o is 
trivial and generates a factor {N~Ni\i — . . . — Ni — M+1). 
The remaining sums are then performed by repeated use 
of the formula 
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For z — Di/Do of absolute value \z\ < 1, we need only 
keep the term i?/(l — z), because in further summa- 
tions the remaining terms generate sub-dominant con- 
tributions of 0(1) or Odzj^). Neglecting these terms 
is justified in our case, because we have assumed that 
the outermost (tunneling) matrix elements are much 
smaller than internal transition elements and by unitar- 
ity, \Di\ < I Do I for all i ^ 0. Keeping only the term 
i?/(l — z) of (|A2|), the structure of the sum always re- 
mains the same, and each summation results in a multi- 
plicative factor (1 — Di/Do)^^. We arrive at the result 
that the sum over paths of length M passing through the 
blocks Do, . . . , Dm is given as 



M 



yN (ii,...,iAf ) 



N D, 



(N-l) 




Con 



1 



Do - A 



C, 



(A3) 
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Eq. (A3) was formulated without aUowing for rep- 
etitions in bfock indices and contains all coupling el- 
ements to lowest order. Loops in block index space 
that stay within the inner blocks can however be in- 
cluded by allowing repeated indices in Eq. ( |A3D . Paths 
with k repetitions of the index combination + 1), 
say, then give rise to contributions containing a factor 
[{Do - D,)-^C'i{Do - Di+i)-^Ci]'' . The number of repe- 
titions can be summed over, which leads to an expression 
as in Eq. (^3) with the replacement 



1 



1 



1 - 



1 



C' 



1 



Do - A+i ' Do-D, 



All types of loops can be included by the corresponding 
replacements, giving rise to a continued fraction structure 



N{ii,...,iM) 
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